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Abstract In this paper we consider two-dimensional, stratified, steady water waves propagating over an 
impermeable flat bed and with a free surface. The motion is assumed to be driven by capillarity (that is, 
surface tension) on the surface and a gravitational force acting on the body of the fluid. We prove the 
existence of global continua of classical solutions that are periodic and traveling. This is accomplished 
by first constructing a 1 -parameter family of laminar flow solutions, 7", then applying bifurcation theory 
methods to obtain local curves of small amplitude solutions branching from T at an eigenvalue of the 
linearized problem. Each solution curve is then continued globally by means of a degree theoretic theorem 
in the spirit of Rabinowitz. Finally, we complement the degree theoretic picture by proving an alternate 
global bifurcation theorem via the analytic continuation method of Dancer. 



1 Introduction 

Initiated by the breakthrough work of Constantin and Strauss (cf. [6]), there has recently been a great deal 
of interest in the theory of traveling water waves with vorticity (see, e.g., [32,33,34,35,36].) Historically, 
of course, most investigations of water waves have focused on the irrotational case, which has several 
extremely convenient mathematical features. While appropriate in a wide variety of scenarios, however, 
the assumption of zero vorticity is far from universally valid. For instance, ocean waves typically have a 
rotational layer near the surface due to the effects of wind shear (cf. Ell ). 

In fact, the wind is known to itself generate steady waves. When it first begins to blow over a body 
of still water, the wind produces a fine, hexagonal meshwork of wavelets whose motion is governed en- 
tirely by capillarity (i.e. surface tension.) If the breeze is maintained, then these waves will grow larger, 
eventually evolving into capillary-gravity waves and finally into gravity waves. Over the course of this 
progression, the waves develop a remarkable robustness: while pure capillary waves will quickly dissi- 
pate in the absence of external forcing, capillary-gravity and gravity waves can propagate long distances 
unaided (cf. EQ). 
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Surface tension's central role in the formation of wind-driven waves argues strongly for the mathe- 
matical consideration of steady, capillary-gravity waves with vorticity, particularly since wind-waves are 
one of the most important classes of rotational water waves. The goal of the present paper, therefore, is to 
develop a global theory for traveling, periodic, water waves with surface tension building upon the local 
theory derived by Wahlen (cf. P2ll3^l331). 

Briefly, our new results come in five parts. The first two of these relate to the phenomenon of double 
(local) bifurcation for capillary-gravity waves. It is well known in the theory of irrotational capillary- 
gravity waves that, when the coefficient of surface tension resides in a distinguished set, the eigenvalue of 
the linearized problem at the point of bifurcation will be of multiplicity two and is otherwise simple. This 
distinctive feature was first observed in the irrotational case by Wilton (cf. 1371 ) and has been the subject 
of numerous papers since. Wahlen observed the presence of double bifurcation points in the rotational case 
and proved the existence of local bifurcation curves in certain regimes (cf. [32]). These results, however, 
are incomplete in the sense that they do not necessarily describe the entire solution set near the bifurcation 
point. The first of our contributions is to describe in detail the local bifurcation diagram obtained in 1321 at 
double bifurcation points. In particular we confirm the existence of curves of mixed solutions conjectured 



by Wahlen (see the discussion following Theorem 1.2 1. Secondly, using the bifurcation theory for analytic 



operators pioneered by Dancer, we extend the local curves to obtain global continua of solutions. The 
application of these techniques to steady rotational waves is in and of itself novel and pays some additional 
dividends. Most notably, it leads to our third result, which is proving that the solution continua are path- 
connected, a fact that was unknown even for the constant density, gravity wave case originally considered 
in [6]. Fourthly, for the case of simple local bifurcation, we also use the more standard degree theoretic 
method to obtain a global bifurcation theorem. Finally, employing the machinery developed by the author 
in ||36ll , all of these results are obtained for fluids with general (stable) stratification, which is entirely new, 
even for the local bifurcation theory. This improvement is quite important from a physical standpoint, 
since stratification phenomenon is both very common in ocean waves (due to salinity and temperature 
variations), and exerts a potentially significant effect on the dynamics of the flow (see, e.g., 1381 . or the 
discussion in 11361 .) 

That said, let us now recall the setup for traveling, stratified, capillary-gravity waves. If we imagine a 
wave on the open ocean, past experience suggests that it may be regular in the following sense. First, it 
is essentially two-dimensional. That is, the motion will be identical along any line that runs parallel to a 
crest. Second, if we regard the wave in a coordinate system moving with some constant speed, it appears 
steady. Finally, the profile is periodic in the direction of motion. 

We now formulate governing equations for waves of this form. Fix a Cartesian coordinate system so 
that the x-axis points in the direction of propagation, and the y-axis is vertical. We assume that the floor 
of the ocean is flat and occurs at y — —d. Let y = rj(x, t) be the free surface at the interface between the 
atmosphere and the fluid. We shall normalize jj by choosing the axes so that the free surface is oscillating 
around the line y — 0. As usual we let u = u(x,y,t) and v = v(x,y,f) denote the horizontal and vertical 
velocities, respectively. Let p = p(x, y, t) > be the density. 

Incompressibility of the fluid is represented mathematically by the requirement that the vector field 
(u, v) be divergence free for all time: 

u x + v y = 0. (1.1) 

Taking the fluid to be inviscid, conservation of mass implies that the density of a fluid particle remains 
constant as it follows the flow. This is expressed by the continuity equation 

p, + up x + vp y = 0. (1.2) 

Next, the conservation of momentum is described by Euler's equations, 

\ pU, + pUU x + pVUy - -P x , 

pV, + pUV X + pVVy = -Py ~gp, 
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y = -d 



Fig. 1.1 The fluid domain D n . The wave is assumed to propagate to the left with speed c. 



where P = P(x, y, t) denotes the pressure and g is the gravitational constant. Here, of course, we assume 
that the only external force acting on the fluid is gravity. 

For capillary waves, we formulate the boundary condition on the free surface as folows. At each point 
we say that the pressure of the fluid deviates from the atmospheric pressure of the air above, Patm, m 
proportion to the (mean) curvature at that point. This results in the dynamic boundary condition 



P = P a 



on y — rj{x, t), 



(1.4) 



(!+>??) 5 

where <x > is the coefficient of surface tension, a material property that we take to be given and constant. 

The corresponding boundary condition for the velocity is motivated by the fact that fluid particles that 
reside on the free surface continue to do so as the flow develops. This observation is manifested in the 
kinematic condition: 

v = 77, + M] x , on y = rj(x, t). (1 .5) 

Since we cannot have any fluid moving normal to the flat bed occurring at y = -d, we require 



0, 



on y 



-d. 



(1.6) 



Note that there is no accompanying condition on u because in the inviscid case we allow for slip, that is, 
nonzero horizontal velocity along solid boundaries. 

We seek traveling periodic wave solutions (u,v,p, P,rf) to ( |l.l[ )-( [L6| l. More precisely, we take this 
to mean that, for fixed c > 0, the solution appears steady in time and periodic in the x-direction when 
observed in a frame that moves with constant speed c to the right. The vector field will thus take the form 
u — u(x - ct,y), v = v(x - ct,y), where each of these is L-periodic in the first coordinate. Likewise for the 
scalar quantities: p(x, y, t) = p(x - ct, y), P — P(x - ct, y), and r\ — jj(x - ct), again with L-periodicity in the 
first coordinate. We therefore take moving coordinates 

(x - ct,y) t-> (x,y), 

which eliminates time dependency from the problem. In the moving frame < |1-1| >-< |1- 3| > become 



u x + v y = 0, 

(U - C)p x + VPy - 0, 

p(u - C)U X + pVUy — —P x , 

p(U - C)V X + pVVy = -Py 



(1.7) 



gp, 
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throughout the fluid domain. Meanwhile, the reformulated boundary conditions are 

v - {u- c)r] x , on y = r]{x), 

v = 0, ony = -d, (1 8) 

P = -Patm - cttj xx (l + i£J , on y = 77C*), 

where m, v,p, P are taken to be functions of x and y, tj is a function of jt, and all of them are L-periodic in 
x. 

In the event that u — c somewhere in the fluid we say that stagnation has occurred. This is a slight 
abuse of terminology, since in our definition we require only that the horizontal velocity vanishes. In order 
to avoid roll-up and other instability phenomena, we shall restrict our attention to the case where u < c 
throughout. 

Recall that we have chosen our axes so that rj oscillates around the line y = 0. Mathematically this 
equates to 

L 

i](x)dx = 0. (1.9) 

In effect, this couples the depth to the problem, so we need not treat d as a free parameter. The trade-off, 
as we shall see, is that this normalization will have significant technical consequences later. 

Observe that by conservation of mass and incompressibility, p is transported and the vector field is 
divergence free. Therefore we may introduce a (relative) pseudo-stream function iff = iff(x,y), defined 
uniquely up to a constant by: 

<A.v = - y/pv, fty = \fpiu - c). 

Here we have the addition of a factor involving p to the typical definition of the stream function for an 
incompressible fluid. This neatly captures the inertial effects of the heterogeneity of the flow (see the 
treatment in 1381 . for example). The particular choice of yfp is merely to simplify algebraically what 
follows. 

It is a straightforward calculation to check that iff is indeed a (relative) stream function in the usual 
sense, i.e. its gradient is orthogonal to the vector field in the moving frame at each point in the fluid 
domain: 

(U - C)iff x + Viffy = 0. 

Moreover, dO) implies that the free surface and flat bed are each level sets of iff. For definiteness we choose 
iff = on the free boundary, so that iff = -po on y — -d, where po is the (relative) pseudo-voumetric mass 
flux: 

yjp(x,y) [u(x,y)~ c]dy. (1.10) 

d 

It is an easy calculation to check that po is well-defined (cf. [ 36 1). Physically, po describes the rate of 
fluid moving through any vertical line in the fluid domain (and with respect to the transformed vector field 
yfp(u-c, v).) 

Since p is transported, it must be constant on the streamlines and hence, we may think of it as a function 
of iff. Abusing notation we may let p : [po, 0] — > R + be given such that 

p(x,y) = p(-iff(x,y)) (1.11) 

throughout the fluid. When there is risk of confusion, we shall refer to the p occurring on the right- 
hand side above as the streamline density function. We shall focus our attention on the case where the 
density is nondecreasing as depth increases. This is entirely reasonable from a physical standpoint. Indeed, 
hydrodynamic stability requires that the depth be monotonically increasing with depth, so that fluids with 
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this type of density distribution are referred to in the literature as stably stratified. The level set -iff — po 
corresponds to the flat bed, and the set where -iff = corresponds to the free surface. Therefore, we require 
that the streamline density function is nonincreasing as a function of -iff. Note that this does not preclude 
the case of constant density. 

By Bernoulli's theorem (for gravity waves) the quantity 

E := P + P -[(u - cf + v 2 ) + gpy, 

is constant along streamlines, that is 

(u - c)E x + vE y - 0. 

This can be verified directly from In the case of an inviscid fluid, E represents the energy of the 

fluid particle at (x,y). The first term on the right-hand side, P, gives the energy due to internal pressure. 
The second and third terms combined describe the kinetic energy, while the last is gravitational potential 
energy. In particular, evaluating the above equation on the level set iff — 0, we find 

|V(A| 2 + 2gp(x, tj(x)) (i](x) + d) = Q, on y = tj(x) (1.12) 

where the constant Q :— 2{E\ n - P a tm + gplijd) gives roughly the energy density along the free surface of 
the fluid. We shall view Q as parameterizing the continuum of solutions. 

Under the assumption that u < c throughout the fluid, and given the fact that E is constant along 
streamlines, there exists a function B : [0, \po\] — > R such that 

^-(x,y) = -Bty(x,y)). (1.13) 
diff 

For want of a better name we shall refer to B as the Bernoulli function corresponding to the flow. Physically 
it describes the variation of specific energy as a function of the streamlines. Define 

B(p):= f B(-s)ds 
Jo 

for po < p < and let B have minimum value S m i n . 

We now briefly outline a few notational conventions. Let D n denote the closure of the fluid domain 

D n := i(x,y) e R 2 : -d < y < r)(x)} . 

For any integer m > 1 and a e (0, 1), we say that a bounded domain D c R 2 is C m+a provided that at each 
point in the boundary, denoted 3D, is locally the graph of a C m+a function. That is, a function with Holder 
continuous derivatives (of exponent a) and up to order m. Furthermore, for fixed m > 1 and a e (0, 1) we 
define the space C^"(D) to consist of those functions / : D — » R with Holder continuous derivatives (of 

exponent a) up to order m and that are L-periodic in the first coordinate. Similarly, we shall take C™ r (Z)) to 
be the space of m-times continuously differentiable functions which are L-periodic in the first coordinate. 
In this paper, we shall seek solutions to the water wave problem (u, v,p, i], Q) in the space 

y := C 2 e + ; (PJ x C^iDJ x C 2 e 7(A;) x C^(R) x R. (1.14) 

Occasionally, for brevity, we shall abuse notation and write (u,v,p,rj) e y, or simply (w,v, rf) e 5f , by 
which we mean (u, v,p, rj, Q) e ^ ', where p is found from the (given) streamline density function using 
< | 1 - 1 1 [ > and Q is calculated via \\.\2\ . 

Finally, in describing the structure of the continua of solutions, we shall make use of the following 
definition. 
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Definition 1.1 A subset "K C 5? is said to admit a global, locally injective, continuous parameterization 
provided that *7C can be described as the image of a locally injective k e C([0, oo); S*). Furthermore, for 
m > 1, we say that K has a locally injective, C m -reparameterization (resp. analytic reparameterization) 
if, for each s« e (0, oo), there exists a continuous t : (—1,1) — » R + such that r(0) = 5* and the map 
5 t-> k(t(s)) is injective and of class C' n (resp. analytic). 

There are several key features distinguishing the sets described above from the standard C m -curves of 
differential geometry. Observe that, while k may be locally injective, Definition 1 1 . 1 1 does not require it to 
be globally injective, i.e., the curve 7C may self-intersect. This scenario is not permissible for C'"-curves 
since in the neighborhood of a point of self-intersection "70 fails to be a one-dimensional submanifold of 
5? . Also, it is important to note that we do not require k' to be non-vanishing. As a consequence, the global 
chart may not be of class C m (resp. analytic), even though at each parameter value we can reparameterize 
to obtain that level of smoothness. 

Our notation established, we may now state our first theorem for stratified capillary-gravity waves. 

Theorem 1.2 (Global Continuum from Simple Bifurcation) Fix a wave speed c > 0, wavelength L > 0, 
relative mass flux po <0 and coefficient of surface tension cr > 0. Fix any a e (0, 1), and let the functions 
8 e C 3+ff ([0, \po\Y) and p e C l+a ([po, 0]; R + ) be given such that the ( |L-B[ ) condition holds (see Definition 



3.5 I. Also we assume the streamline density function p is nonincreasing. There exists a set Z\ c (0, oo) 



such that the following is true. E\ contains a half-line, (cr c , oo), and, if a e E\, then we have: 

(a) There exists a connected set C C 5? of solutions (u,v,p,rj, Q) of the traveling, stratified, water wave 
problem with surface tension ( |l.l| l-( fL6} . The setC, moreover, has the following properties. 

• C contains a laminar flow (with aflat surface jj = and all streamlines parallel to the bed); denote it 
byiU .() ). 

• Either (i) C contains more than one laminar flow, or it is unbounded in the sense that, along some 
sequence {(u n ,v n ,p n ,T] n )} c C, we have 

(ii) max u n \ c or (iii) min u„ J, -oo. 

(b) Furthermore, there exists a path-connected subset 7C c C such that 

• 7C admits a global, locally injective, continuous parameterization with a locally injective C l reparam- 
eterization. 

• Either 7C is a closed loop, or it unbounded in the sense that alternative (a)(ii) or (a)(iii) must occur 
along "70. 

(c) All flows (u, v,p, T]) e C are regular in the sense that: 

• u, v, p and rj each have period L in x; 

• u, p and r] are symmetric, v antisymmetric across the line x — 0. 



Remark 1.3 The Local Bifurcation Condition ( |L-B| l essentially states that there exists an eigenvalue of the 
linearized problem with lowest mode, and is both necessary and sufficient for our result to hold. Still, in 
Lemma 3.6 we derive an explicit sufficient condition that implies ( |L-B[ ): 



(gp(0) + a)pl > 



f 

•J Pa 



4tt 2 



(2B(p) - 2B min + 2e ) 3/2 



+ (p - po) 2 ((2B(p) - 2B min + 2e ) 1/2 + gp'ip)) }dp, 



(1.15) 
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where €q is defined by 



3/2 
fc 



max {2g\\p'\\ L ~pleM, (2g\\p'\\ L ~? , (4||p'[| L ») 3 , Sg\p \p(O)} . (1.16) 



This is most certainly not optimal. Indeed, obtaining a complete characterization of the set of data {po,p,f3) 
for which ( |L-B[ ) holds is an important and largely open problem, even in the special case when p = 1 and 
cr = 0. It is worth mentioning that ( |1.15| l can be satisfied by assuming \po\ and ||p'||i» are sufficiently small, 
or cr is large. From this it is easy to see that the set E :— {cr > : ( |L-B| l holds} contains a neighborhood of 
+oo. 

We note also that much is known about the structure of the continuum C near the point of bifurcation 
([/„, (?*)• There exists a neighborhood H of ([/,, Q*) in 5? in which CC\H - 'Knti comprises all solutions 
to ( |l.l| l-( [L6| , and (U„,Q t ) is the only laminar solution in H. Moreover, all nonlaminar flows in 1/ n C 
exhibit the following monotonicity properties: 

• within each period the wave profile i] has a single crest and trough; say the crest occurs at x — 0; 

• a water particle located at (x, y) with < x < L/2 and y > -d has positive vertical velocity v > 0. 

For the proof of these results, see Theorem |3.1| 

Finally, let us note that in certain special cases it can be shown that2i is a countable set, see Proposition 

□ 



As we shall see, Theorem 1 .2 is obtained by means of a bifurcation argument with the family of laminar 
flows taking on the role of the trivial solutions. When cr e E\ (for example, if cr is sufficiently large) we 
will show that local bifurcation occurs at a simple (generalized) eigenvalue. 

On the other hand, for cr £ Z\, the eigenvalue will be not be simple. In order to discuss this scenario 
more carefully, we must first make a few definitions. Let D be a smooth domain in R 2 and let Banach 
spaces X,Y a CT£ a (D) be given for some m > 1, a e (0, 1). For a smooth map T : R X X — > Y with 
0) = 0, the generic bifurcation problem is to find all nontrivial solutions (A, w) of 

r(/l,w) = 0, (1.17) 

in a neighborhood of a given solution (A„, 0). Suppose that the null-space of T w {A t ,Q) is two-dimensional, 
spanned by {(/> \ , <p2 }, and <p, = <pi(x, y) is of minimal period L / n, in the first coordinate for i = 1,2, where n \ 
and «2 are distinct natural numbers. Assume further that there exists a (continuously) parameterized curve 
of solutions, 

Cioc := {(A(s),w(s)) : s e (-e, f))clx X, 
to ( |1.17[ l with (A(0), w(0)) = (A„ 0). We shall refer to these solutions as pure provided that 

w(s) = s<pi + o(s) in X, 

for either i = 1 or i = 2 but not both. That is, pure solutions are those that locally lie near the interior 
of the span of either cp\ or 02. Conversely, the elements of C\ oc are called mixed solutions if there exists 
fi,£> e R x such that 

w(s) = s (£i0i + £ 2 02) + o(s) in X. 

Clearly, mixed solutions locally lie near the interior of the span of <p\ and <p2- 

The relevance of these concepts to the case at hand stems from the fact that, for cr e R + \ E\, local 



bifurcation occurs at eigenvalues of multiplicity two (see Corollary 3.14 ) That is, the null space of the 
linearized operator is spanned by two functions, {<p\, fc], of minimal period L/rii,L/ri2 in x, respectively. 
This feature of the eigenvalue problem is well-known for homogeneous capillary-gravity waves and has 
been very well analyzed in the irrotational case (cf. [15 25,30 1 for analytic results and, e.g. [26 28 29 1 
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and the references contained therein for some numerical results.) For general vorticity, this scenario was 
analyzed by Wahlen [32|. Briefly, his argument was the following. So long as there is no resonance, 
i.e. «2 is not an integer multiple of n\, then the dimension of the null space can be reduced to one by 
considering the restriction to the subspaces of Ljn\- or L/n2-periodic solutions. Applying the Crandall- 
Rabinowitz arguments developed for the simple eigenvalue case, one can prove the existence of two C 1 - 
curves bifurcating from the family of laminar flows, say Cij c an d C2,ioc- These consist of small amplitude 
solutions of minimal period L/tii and hjni, respectively. In fact, Cy c are both curves of pure solutions in 
the sense outlined above. 

The price payed for the relative ease of this approach is that, in restricting the solutions space, we 
are likely missing the mixed solutions. To see why this may be the case, note that a linear combination 
£i$>i + ^202 ma Y have minimal period L. Such solutions would therefore not be included in either Ci i oc or 

C2,loc- 

Worse still, when resonance does occur (i.e. niln\ £ N), this type of argument offers even less infor- 
mation. Since restriction to the L/rii -periodic solutions no longer renders the null space one-dimensional 
(as L/«2-periodic solutions are now also L/nj -periodic), the only conclusion we can draw via Crandall- 
Rabinowitz is that there is at least one local curve of pure L/«2-minimal periodic solutions. In the language 
we have been using thus far, this means that we can prove the existence of C2,ioc> but not Cijoc- This, of 
course, says nothing about the existence (or nonexistence) of curves of mixed solutions. 

In the present paper we devote some effort to advancing the study of the bifurcation near double eigen- 
values. Rather than using Crandall-Rabinowitz, we instead carry out our analysis via the full Lyapunov- 
Schmidt procedure. Owing to a great deal of degeneracy in the problem, a complete account is, unfor- 
tunately, beyond the scope of this paper. Nonetheless, we are able to obtain some detailed results on the 
structure of the full solution set. 

Theorem 1.4 (Local Continua from Double Bifurcation) Let the wave speed c > 0, the relative mass flux 
po < 0, surface tension constant o~, streamline density function p e C 2+a ([0, \po\]) and Bernoulli function 
B e C l+a ([0, \po\Y), < a < 1 be given. Suppose further that po, p, B and cr collectively satisfy ( |L-B[ l. 



In contrast to Theorem 1 .2 suppose cr <£ X,\ but non-degeneracy conditions (HTb and (|4.2[l hold. Then 



there exists four C -curves, (Cy c}; = i> of small- amplitude, traveling wave solutions of ( |l.l| l-( fL6} in the 
space 5? with the following properties. 

• Each (u, v,p, j], Q) e I J j Cij c has speed c, relative mass flux po and satisfies u < c throughout the fluid. 

• Fori- 1,2, each (u,v,p,7],Q) e C^ioc, i - 1,2, has minimal period 2n / «; in x. 

• For i — 1 , 2, 3, 4, each of the curves {Cijoc} contain precisely one laminar flow, while for each nontrivial 
solution (u, v,p, r\, Q) e C u ] oc , with i = 1,2, we have 

(i) the functions u and rj are symmetric around the line x = while v is antisymmetric; 

(ii) the function tj has precisely one maximum (crest) and one minimum (trough) per minimal period; 
and 

(iii) the wave profile is strictly monotone between crest and trough. 

Remark 1.5 Referring back to the preceding discussion, the first two curves, Ci j oc and C2,i c are the pure 
solutions, while the remaining two are the mixed solutions. 

The nondegeneracy conditions ( |4.1| l-( |4~2"| ) ensure that, upon applying the Lyapunov-Schmidt reduc- 
tion, the solutions of the bifurcation equation are described totally by the expansion to cubic order. It is 
possible to omit condition ( |4.2[ ) but, in this case, rather than having the complete solution set near ([/*, Q„), 
we can only guarantee the existence of TCi . Condition ( |4.1| l can be similarly relaxed (see the remarks in 
section^. More detailed information on the local structure of the bifurcation diagram is outlined in Lemma 
|4!21-Lemma 14771 and Theorem gTS] □ 

Once the local bifurcation picture has been resolved, we are able to continue the curves {C,j oc }^ =1 to 
produce global continua of solutions. All told, this gives the last of our main results: 
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Theorem 1.6 (Global Continua from Double Bifurcation) Assume the hypotheses of Theorem 1 .4 There 
exist global continua of solutions extending {Ci,\ 0C } A i=x in the sense that the following statements hold. 

(a) There exists a path-connected set 7C c 5? of solutions (u, v,p, rj, Q) to \\.\ | )— ( p~6] >. The set "70 has the 
following structural properties. 

• 7C contains a laminar flow, ([/*, Q„). 

• There exists four path-connected sets {'Ki}*^ such thafK — Ui^C' an d (U*, Q„) e Hi^C'- in fact, there 
exists a neighborhood U of(U*, Q*) in J? such that 

Qioc = % nl/, fori = l, 4. 

• Each *Ki admits a global, locally injective, continuous parameterization with a locally injective C l - 
reparameterization. 

• For each i = 1, . . . , 4, we have the following alternatives: Either (i) % is a closed loop or, along some 
sequence {(u n ,v n ,p„,n n , Q n )} c %, 

(ii) max u n \ c or (iii) min u n J, -co. 

(b) All flows (u, v,p, rj) e 7C are regular in the sense that: 

• u, v, p and r\ each have period L in x; 

• u, p and rj are symmetric, v antisymmetric across the line x — 0. 

Steady capillary-gravity waves, as we remarked earlier, are an important, classical object of study in 
fluid mechanics and geophysical fluid dynamics. Unsurprisingly, therefore, they have been the subject of 
a truly voluminous quantity of research, far more than we can hope to completely summarize here. The 
vast majority of the work on steady, capillary-gravity waves considers the irrotational regime. There are 
physical justifications for this choice, but perhaps the most compelling reason is mathematical: when there 
is no vorticity, tff is a harmonic function. All that remains is to determine the free surface n. As the vector 
field has zero curl, moreover, we have the existence of a velocity potential cp, which enables the use of 
conformal maps and other powerful techniques from complex analysis. For instance, there are a number 
of equivalent reformulations of (jTTTJ — (|TT7j as a scalar equation for 6, the slope of i] (cf. [ 20 24 j for the 
classical case where <x = 0, and, e.g. [25 1 for the generalization to <x > 0.) 

Of course, with the introduction of vorticity, the usefulness of these tools is greatly diminished. For 
that reason, the structure of our arguments in this paper follows the general program for steady, rotational 
waves established in [6 33 32 36] . We begin, in section|2] by reformulating the problem in terms of scalar 
quantities and in a semi-Lagrangian coordinate system. As a consequence, the transformed domain will 
be fixed, and the system of equations ( |l.l[ )-( fL6| is replaced by a non-linear, integro-differential, boundary 
value problem for a scalar function, h, that describes the height above the flat bed. The interior equation 
satisfied by h has the following form: 

£(D 2 h, Dh) + I(Dh, h) = 0, 

where D 2 h denotes the Hessian of h, Dh the gradient, X is an elliptic operator and I is an integro- 
differential operator. As in 11361 . we will show that I totally accounts for the presence of stratification 
(indeed, if p is constant, then I = 0). On the other hand, the technical challenges due to surface tension 
present themselves, naturally, on the boundary. The reformulation of the dynamic condition (in particular, 
the mean curvature term in ( |1.4[ |) results in a non-degenerate Venttsel-type condition on the image of the 
free surface. 

In section [3] we prove the existence of a 1 -parameter family, 7", of laminar flows. Since these flows 
are, by definition, functions of the streamlines, the integro-differential equation reduces to an ODE with 
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operator coefficients (the operator coming from I). We proceed to show that, if cr e E\ — e.g., it is suffi- 
ciently large — then there exists a point on T at which the problem has a simple (generalized) eigenvalue. 
Using the theory of Crandall-Rabinowitz, we prove that there is a branch of small- amplitude, non-laminar 
solutions bifurcating from T at that point. This is, essentially, a generalization of [32| to the stratified case. 

Next, in section [4] we consider the case when cr £ E\. As mentioned above, this means that the 
bifurcation from T is not simple. In fact, we show that, for all but at most one value of cr, the bifurcation 
is precisely double, i.e. the null space of the linearized operator has algebraic multiplicity two. Using the 
classical method of Lyapunov-Schmidt, we give several results on the bifurcation diagram for such cr. This 
is done by explicit computation of the coefficients in the bifurcation equation. In particular, we show that, 
under the non-degeneracy condition, there are four C 1 -curves of small-amplitude, non-laminar solutions 
emanating from the bifurcation point, as described in Theorem 1.4 Let us note that this result, though 
preliminary in the sense that it requires some additional non-degeneracy assumption, is new even in the 
case of constant density. 

The remainder of the paper is devoted to the global analysis. In section [5] we use a degree theoretic 
argument to prove an alternative theorem in the style of Rabinowitz [27]. This allows us to continue 
globally the curve found via simple local bifurcation in section [3] Of course, this involves first proving 
some strong compactness properties about the reformulated problem — a task made difficult both by the 
nonlocal operator term in the interior equation, and the nonstandard, second-order boundary conditions. 

To continue the local curves obtained in the case of double bifurcation, we employ a global bifurcation 
theory due to Dancer (cf. [5,9 10, 1 1 J ) that exploits the analyticity of the integro-differential operator in 
the reformulated problem. As we make clear, this method enjoys several advantages with respect to the 
more standard degree theoretic approach of section [5] First, it enables us to prove an alternative theorem 
that continues each of the local curves found in section |4] Since the crossing at the bifurcation point is no 
longer odd, this would be non-trivial with a degree theoretic argument. Second, we are able to prove higher 
regularity of the global continuum. In fact, we show that the continuum "70 of solutions to the reformulated 
problem is a curve with a locally injective, analytic reparameterization. 

In section [7] we derive uniform bounds on the C 3+ff -norm along the continuum C in terms of the 
first-order derivatives, the mean curvature on the surface and the energy Q. In section [8] these bounds, 
along with the alternative theorems of sections [5] and [6] are combined to produce Theorem |1.2| and The- 
orem |1.6| Finally, section [9] analyzes some important special cases: homogeneous, irrotational waves and 
homogeneous rotational waves. 



2 Reformulation of the problem 

In this section we set ourselves to the task of finding an equivalent formulation of ( |1.7| l-( fL"8j ) that elimi- 
nates the free boundary problem. We shall accomplish this by making a change of variables to transform 
the fluid domain in a fixed domain, whose closure we shall denote R. As usual for incompressible 
problems, this will be done by means of the pseudo-stream function i/c and will come at the cost of some 
additional nonlinearity in the governing equation. 

To see why this is a natural choice of coordinates, recall that we have ft = on the free surface, 
iff = -p on the flat bed, where po is the relative pseudo-mass flux and ip is defined uniquely by the 
requirement that 

<A.v = - yfpv, if/y = y[p(u - c). (2.1) 
In light of \2.\ \ and (JO) -(TO), the governing equations inside the fluid become 

=~ p * in 5; (2.2) 



Water waves with surface tension 



11 



whereas the boundary conditions (|1.8[) are 



' \j/ x = -i/fytlx, on y = ?7(x), 

P = P a tm ~ o-7] xx (l + Til) 3/2 : , on y = rj(x), (2.3) 
(A v = 0, on y - -d. 

Also, by Bernoulli's theorem, the quantity 

E = P+^((u-c) 2 + v 2 )+gpy, (2.4) 

is constant along streamlines. In particular, then, evaluating the above equation on the free surface we find 

|V<A| 2 + 2gp(x, rj(x)) (r](x) +d) = Q, ony = rfa), (2.5) 

where the constant Q :— 2(E\ n - P atm + gp\rjd). Note that this Q gives roughly the energy density along the 
free surface of the fluid. Critically, however, we shall see that as a consequence of our normalization of t], 
d is not a parameter for the problem. On the contrary, in all but the most trivial cases, d varies along C. In 
our analysis, therefore, we shall instead be viewing Q as parameterizing the continuum. 

With that in mind, we begin by deriving a scalar PDE satisfied by if/. In the case where the fluid is 
homogeneous, it is not hard to check that 

-Atff = a>, (2.6) 

where a> :- v x — u y is the vorticity. Moreover, again exploiting the homogeneity, one can readily verify 
that oj is transported, i.e. there exists some vorticity function, y : [0, \po\\ — » R, such that a> - y(xji). 
It is therefore possible to view ( |2.6| l as a semilinear elliptic equation, where y is prescribed. This is the 
approach taken, for instance, by Constantin and Strauss ((6,71) in the (homogeneous) gravity wave case, 
and Warden ([33 , 32 1) for (homogeneous) capillary and (homogeneous) capillary-gravity waves. 

In the presence of stratification, however, we cannot generally expect the vorticity to be a function of 
if/. Instead, we must look back to the Bernoulli theorem to find an appropriate quantity conserved along 
streamlines. In particular, one can prove that 

dE ^ ^ dp 
dif/ ^ dif/ 

Recall that in the previous section we introduced the Bernoulli function /3 and streamline density function 
p. Rewriting the above expression we arrive at an equation of roughly the same form as (|2.6]l: 



-PM=Ai/,-gyp'(-M. (2.7) 

This is known as Yih's equation or the Yih-Long equation (see, e.g., l38l ). In this paper, it shall serve as 
our governing equation for iff. Were it not for the presence of y on the right-hand side, this would reduce 
to the homogeneous case — a fact that agrees with our intuitive notion that density stratification should 
reintroduce the depth into the problem as a serious consideration. One can easily check, moreover, that, if 
p = 1 , then f3 and y coincide exactly. 

With ( |2.7| l in hand, we now make a change of variables to eliminate the free boundary. The new 
coordinates we will denote {q, p) where 

q = x, p = -i/f(x,y). 

This scheme is sometimes referred to as semi-Lagrangian coordinates, in recognition of the fact that we are 
working, in some sense, halfway between Lagrangian (or streamline) coordinates and the usual Eulerian 
system (cf. 13T))- 
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By means of a scaling argument, we may take L :- 2n. Then, under the transformation 

(x,y) h> (q,p), 
the closed fluid domain D n is mapped to the rectangle 

R := {(9, p) e R 2 : < q < 2n, p < p < 0}. 
Given this, it will be convenient to put 

T := {(?,/>) e R : p = 0), 5 := e R : p = p }. 

Note that, in light of ( |l.ll| l and ( |1.13| l, we have that 

P=ft-p), P=P<J>). 

Moreover, the assumption that the streamline density function is nonincreasing becomes 

P„ < 0. (2.8) 

Next, following the ideas of Dubreil-Jacotin (cf. 11211131 '). define 

h(q,p):=y + d (2.9) 

which gives the height above the flat bottom on the streamline corresponding to p and at x = q. We 
calculate: 

1 K 

xlf y = --, <A. V = (2.10) 
h p h p 

Note that this implies h p > 0, because we have stipulated that u < c throughout the fluid. The change of 
variables then gives 

1 K 
u = c-——, v = -— 1-. (2.11) 
yjph p y/ph p 



Performing the Dubreil-Jacotin transformation on ( |2.7[ l, we find that, in the interior of the fluid, 

(1 + h 2 q )h pp + h qq h 2 p - 2h q h p hp q - g{h - d(h)) P ph\ = -h)fi(-p), 



where we have used ( |2.9| l to write y — h - d. Recall, however, that we have normalized 77 so that it has 
mean zero. Taking the mean of ( |2.9[ ) along T, we obtain 

pin 

d = d(h) = \ h(q,0)dq. (2.12) 
Jo 

That is, the average depth d must be viewed as a linear operator acting on h. Namely, it is the average value 
of h over T. Where there is no risk of confusion, we shall suppress this dependency and simply write d. As 
we shall see, the addition of this integral term to the governing equations will be the single most significant 
departure from the homogeneous case, both technically and qualitatively. In recognition of this fact, we 
shall refer to the mapping h t-> -g(h - d(h))h 3 p p p as the stratification operator. See ||36ll for details. 

On the bottom, we have by construction that h = 0. The more interesting boundary condition is found 
on the top of R, which corresponds to the free surface in the original formulation. There, we have from 
( [23] ) that 

1 + h\ + h 2 p (2o- K [h] + 2gph -0 = 0, p = 0, 
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where the curvature, k, is defined by 



hqq 



3/2 ' 



This calculation is carried out in l33l . for example. 

Combining these observations together we see that the completely reformulated problem is the fol- 
lowing. For a given Bernoulli function B, find (h, Q) e C^"(R) X R with h p > and satisfying the height 
equation 

' (1 + htyhpp + h qq h 2 p - 2h q h p h pq 

-g(h - d(h)) Pp h 3 p = -hfa-p), p <p<0, 

(2.13) 

1 + h\ + h 2 p (2<TK[h] + Igph -Q) = 0,p = 0, 
h = 0, P=Po- 

Lemma 2.1 (Equivalency) Problem ( |2.13| l is equivalent to problem 



Proof This is proven for cr - in [36|. A straightforward modification of that argument — which is itself 
a simple generalization of a similar result in |6) — proves the lemma. For brevity, we omit the details here. 
The reader is also referred to the discussion in Remark |8T1 □ 



3 Local bifurcation at simple eigenvalues 

3.1 Overview 

The goal of this section is to prove the existence of local curves of small amplitude solutions to (|1.7[ )-( [l.8| > 



In continuation from the previous section, we shall make use of the equivalent formulation p. 13 i to prove 



the existence of a 1 -parameter family of laminar flow solutions, T . These will serve as the trivial solution 
curve for the standard Crandall-Rabinowitz machinery. Linearizing the problem around the laminar flows, 
we are thereby able to prove the existence of a local curve of non-laminar solutions bifurcating from T at 
a simple eigenvalue under certain size conditions. The analysis culminates in the following result. 

Theorem 3.1 (Local Bifurcation at Simple Eigenvalues) Fix a e (0, 1) and let the wave speed c > 0, the 
relative mass flux po < 0, surface tension constant cr, streamline density function p e C 2+a ([po, 0]; R + ) 
and Bernoulli function B e C l+a ([0, \po\J) be given. Suppose further that po, p, B and cr collectively satisfy 
( |L-B| >. IfcreZ\ (as defined in ( |3.6| l), then the following holds. 

There exists a C ] -curve Ci oc c 5? of small amplitude, traveling wave solutions (u, v,p, rf) of ( |l.l[ l-( [L6| 
with period 2n, speed c and relative mass flux po, satisfying u < c throughout the fluid. The curve Ci oc has 
the following properties. 

(i) Cioc contains precisely one laminar flow, denote it (Q*, U*)- 

(ii) Cioc \ KG*, Ut)} comprises all nonlaminar solutions in a sufficiently small neighborhood of(Q t , [/*) in 

Rxc 2 P : r a (T^). 

(iii) For each nontrivial solution (u, v,p, rf) e Ci oc , 

• the functions u, p and t] are symmetric around the line x — while v is antisymmetric; 

• the function rj has precisely one maximum (crest) and one minimum (trough) per period; 

• the wave profile is strictly monotone between crest and trough. 
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3.2 Eigenvalue problem 



We begin by establishing the existence of a 1 -parameter family of laminar flows. By laminar flow we 
mean that the streamlines (which, we recall, include the free surface) are parallel to the bed. Waves of this 
type will accordingly have the ansatz h = H(p), which reduces the height equation ( |2.13| l to the following 
(operator coefficient) ODE boundary value problem: 



H pp - gHj(H - d{H))p P = -H;/3(-p), p <p<0, 

l+H 2 p (2gpH-Q)=0, p = Q, (3.1) 

H = 0, p = po. 



The basic existence theory for laminar solutions comes in the form of the following lemma. 

Lemma 3.2 (Laminar Flow) Suppose that the streamline density function p satisfies p p < 0. Then there 
exists a 1 -parameter family of solutions (H(-;A), Q(A)) to the laminar flow equation with H p > 0, 
where < -2B min < A < Q. 

Proof Because the curvature of the laminar flows is zero, ( |3.1| ) is identical to the problem satisfied by 
laminar stratified gravity waves in the absence of surface tension. The lemma therefore follows from 
Lemma 3.2 of [361. □ 



Remark 3.3 For later reference, we note it was also shown in l36l that 

H(0) = T~A' H n (0) = A ' m - 
2gp(0) 

In contrast to the theory for homogeneous gravity and capillary-gravity waves, there is no explicit 
formula for these laminar flows. However, H p satisfies the following implicit equation 

H p (p; A) = 1 -, po < p < 0, 

where 

G(p; A) : = 2B(p) + 2 f g (H(r; A) - d(H(r; A))) p p (r)dr, P o<P<0. 
Jp 

It will occasionally be convenient to make use of the shorthand 

a(--A) := H p (-;A)-\ Y(-;A) := //(•; A) - d{H). 



Finally, we shall denote H := dH/dA, Y := dY/dA, G := dG/dA and Q := dQ/dA. In fact, we have the 
following useful identity 

np) - \ £ iT7GWr> dr = U° ( J + ^ air) ~ 3dr > p ^ p ^°- (3 ' 2) 

This is found by differentiating the relation Y„ — (A + G) 1 ^ 2 in A, followed by integrating in p. □ 
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We now linearize the full height equation ( |2.13| l along the curve of laminar solutions, which we shall 
denote T . As in ll36ll . when there is heterogeneity in the fluid, it will be necessary to ensure that A is 
bounded strictly away from -2B mm . Recall that we have defined 



Note that having e^ 2 greater than the second quantity on the right-hand side above implies, for A > 



^ /2 : =max{2£||p'|L^ e l P ol (2g\\p'\\ L ^\ (4||p'|| L ~) 3 , 8*|po|p(0)}. 

ving e 

-25 m in + eo, 

a + gp P > 4' 2 ~ g\\Pph- > £o/2 > 
if p p S 0. The same holds true, of course, when p p = 0, since in this case the left-hand side reduces to a, 
which is strictly positive when A > -2B mm . It was shown in l36l that, when e^ 2 > 2g\\p p \\ L ™p 2 ) e [Pa<[ , then 
-1/2 < G(-; A) < 0. In particular, this implies that, a+gp p > 0, and 1 +G > 1/2 > 0. Both of these facts are 
critical for our analysis and so from this section onward we shall always be assuming that A > -2B„^ a + eg- 
The remaining two quantities on the right-hand side above are, in fact, not strictly necessary. In [ 36] they 
were used to avoid a certain pathology that might occur if the point of bifurcation A, occurred to the right 
of a particular value Aq. That will not be an issue for the present paper. Nonetheless, for simplicity of 
exposition, we keep the same definition of eo- 

Now, fixing any such A, we calculate that the linearized problem about //(•; A) is: 

V + m qq H 2 - gp p (3m p (H - d{H))H 2 p + (m - d(m))H 3 p ), 

= -3H 2 p mpP(-p) p <p<0, 
m = 0, p = po, 

, gpm - crm qq - A 3/2 m p = 0, p = 0. 

Since we seek solutions that are 2^-periodic and even in q, it is natural to search for m with the ansatz 
m(q, p) — M(p) cos (nq), for some n > 0. Inserting this into ( |3 .3 1 > we see immediately that M must satisfy 
the following: 

- n 2 MHl - gp p (3Mp(H - d(H))H 2 p + (M - <5„ M(0))/^) 
= -3H 2 p MpB(-p), p Q <p<0, 
M = 0, p = po, K ' 

v gpM + n 2 crM - A 3/2 Mp = p = 0, 

where 6,,o = 1 if n = and vanishes otherwise. This is a simple consequence of the fact that 

pin r*2n 

d(m) = -r m(q, 0)dq = M(0) -r cos (nq)dq = 6„ M(0). 
Jo Jo 

Using the equation satisfied by H, ( |3. 1 1 >, we may put ( |3.4| l into self-adjoint form: 

f \HfMp\p = (n 2 H' p l + g Pp )M - gp p 6 n0 M(0), Po <p<0, 

\M = 0, p = po, (3.5) 

[ A 3/2 M p = (n 2 cr + gp)M, p = 0. 

Solving the above problem for any value of n will produce a 27r-periodic solution to the linearized problem, 
but, as we are primarily interested in solutions of minimal period 2k, we shall focus on the case where 
n = l. 

With that in mind, we make the following definition. 

Z\ := {<r e R + : there exists a unique A > -2B min + e such that 
problem \3.5\ has a solution (M, n) if and only if n = 1}, 
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E% '.— {cr e R + : there exists a A > -2B mm + £q such that 

problem (|3.5|l has precisely two linearly independent 



(3.7) 

nontrivial solutions {M\,ti\) and (M 2 ,n 2 ), 
where n\ — \,n-i> 2}. 

and 

E?> : = {cr e R + : there exists a A > -2B m ; n + €q such that 

problem ( |3.5| l has at least two linearly independent solutions (3.8) 
(Mi,«i) and (M 2 ,n 2 ) where n\ = 1, n 2 = 0} 

For fixed po, p and p\ we shall prove that the set E\ consists of those values of cr for which simple 
bifurcation occurs at some point along T . The curve of nonlaminar solutions will, locally, be of the form 
(Q(A), H(p; Q{A)) + eM(p) cos q + 0(e 2 )), and hence of minimal period 2n. On the other hand, if cr lies 
in the set E%, then we expect double bifurcation to occur and there will be several curves of 2^-periodic 
solutions that, nevertheless, have a minimal period strictly less than 2n. The last set, E3, consists of those 
values of cr where there is a zero eigenvalue of the linearized problem. We shall prove that E3 is either a 



singleton or the empty set (see Lemma 3.10 1. Furthermore, in Corollary 3.14 we show that every cr for 
which ( |L-B| > holds is in precisely one of S\, E 2 and Ej . We strongly suspect that E\ is generic. In fact, in 
section 19] we provide an elementary proof of the (previously known) result that, for irrotational, constant 
density waves, E 2 is countable. 

This section is concerned with the simple bifurcation case, so our first task is to produce sufficient 
conditions under which cr e E\. In order to apply the linear theory developed in 061 . we introduce a 
modified density function, p, defined by 

cr 

p :=p+ -. 

8 

It follows that, for each A > -2S m ; n + eo, (H(-; A), Q(A)) is a laminar solution with streamline density 
function p, where 



Q(A):=2crH(0;A) + Q(A) = cr 



'(cr- l gp(0) + l)Q(A)-r 



8P(0) 



gp(0)Q(A) - crA 

gp(0) ■ ^ } 



It is important to note that Q inherits all of the important properties of Q. These we summarize in the 
following lemma. 



Lemma 3.4 Let A > -2S„„„ + eo be given and let Q be defined as in ( |3.9| l. Then Q is strictly convex as a 
function of A with a unique minimum occurring at Aq. Moreover, if the unique minimum of Q occurs at A§, 
then Aq < Aq. 

Proof We know by Corollary 3.7 of [36] that Q is a strictly convex function of A. Differentiatin g (|3.9| > 
twice in A confirms that the same is true of Q. To see Q attains its minimum, note that evaluating ( |3.9| > at 
A — Aq, we have dQ/dA < 0. But, since dQ/dA /* 1 as A — » 00, while gp > cr, taking A sufficiently large 
we must have dQ/dA > 0. Hence, there exists some /lo > ^0 where dQ/dA — 0. □ 

Now, since we are searching for a solution to ( |3.5[ > with n 2 = 1, we can simply replace p with p to 
rewrite the problem as something reminiscent of the stratified gravity wave case. In particular, we see that 
it is equivalent to finding a solution M of 

(H p {H- p 3 M p } p - gp p H p M = M, po < p < 

\M = 0, p = p Q (3.10) 

{A V2 M p =gpM, p = 0. 
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We have therefore reduced the problem to that considered in 061 . Recalling the methodology of that 
paper, let us now consider the minimization problem: 

H = p,(A) = mf <R((f>;A) (3.11) 

tpeS 

where iS := {0 e H ((po, 0)) : £ 0, (f>(po) = 0} and the Rayleigh quotient *R is given by 

-gp(0)m 2 + f° H p (p; A)~ 3 cf> p (p) 2 dp 

<R(<P;A) := —5 ^ . (3.12) 

j po cf>(p) 2 (H p (p;Ar l +gp P (p))dp 

By the standard variational theory for Sturm-Liouville problems, the existence of a nontrivial solution to 
( |3.10[ ) for some A t > -2B mm + eo is equivalent to the statement p(A t ) = — 1. Given that fact, we make the 
following definition. 

Definition 3.5 We say the pseudo- volumetric mass flux po, coefficient of surface tension <x, streamline 
density function p and Bernoulli function B collectively satisfy the capillary-gravity local bifurcation 
condition provided that 

there exists A > -2B m [ n + £q and Mi solving the linearized problem ( |3.10[ >. (L-B) 

Equivalently, 



inf 12(A) < -1. 

,i>-2B min +6o 



This last statement can also be written 



(-gp(0)cf>(0) 2 ~ o-cf>(0) 2 + £ H p (p; AT 3 cf> p (p) 2 dp) 
inf inn 5 ^ \ < -1. (3.13) 

A>-2B mm+ e fieS ^ ^2 ^ Hp(p . A) -l + gpp(p) } dp j 

where S :- {cf> £ H i ((p Q , 0)) : i 0, (p(po) = 0} and //(•; A) is the laminar flow solution given by Lemma 
( |L-B[ > is both necessary and sufficient, but clearly suffers from being highly non-explicit. However. 



as was alluded to in the remarks following Theorem 1 .2 if we are willing to sacrifice necessity, there is 



an entirely explicit sufficient condition under which the capillary-gravity local bifurcation holds, namely 

Lemma 3.6 (Sufficiency of Size Condition) If po, p, B and a satisfy the size condition \\.\5\ , then they 
satisfy ( |L-B[ ). 



Proof This lemma follows from Lemma 3.5 of |36| with p in place of p. □ 

We also remark that from ( |3.13| l it is evident that, if ( |L-B| l holds for (po,cri,p,B), then it hold for 
(po,o- 2 ,p,B) whenever cr 2 > o~\. In particular, this implies that if the Local Bifurcation Condition is 
satisfied for the problem without surface tension (i.e., with cr = 0), then it is satisfied for all values of 
<x > as well. 

For future reference we now quote another result from 11361 . 
Lemma 3.7 (Monotonicity) p is a strictly increasing function of A when p(A) < 0. 
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This result implies, for instance, that any A > -2B m i n + €q satisfying p(A) — - 1 is necessarily unique. The 
next lemma gives the existence of such a A. 

Lemma 3.8 (Eigenvalue Problem) Assuming ( |L-B| >, there exists A t > -2B mm + e§anda nontrivial solution 
M = M(p) to (|3~3)l with n = 1. 



Proof In the previous paragraphs we have shown that the producing such an M is equivalent to finding 
A t such that p(A t ) = — 1. Condition ( |L-B[ ) ensures that p(A) < -1, for some A > -2B mm + €q, so it only 
remains to show that p(A) > - 1 for A sufficiently large. With that in mind, let A be given such that 

A > -2B mm + (g\\p p \\zr + V#P(0) + of . 
Recalling that G(-; A) > 2B mm , we find that, for any such A, 



H- { > H p l + gp p = (A + G) 1/2 + gp p >(A + 2B min ) 1/2 - g\\p p \\ L <» > yjgp(0) + a. 
Now let e 5 be given and take A as above. Then 

(h- 3 ^), + (H~ l + gpp)^ 2 ) dp > V^P(O) + o- ((gp(0) + <r)<f>l + 2 ) dp 

J Pa Jpo 

r° 

>2(gp(0) + o-) <p<p p dp = {gp(0) + o-)m 2 

•Jpo 



From the last inequality we conclude < R((/>) > - 1 . Since the selection of above was arbitrary, this proves 
p(A) > -1, which completes the lemma. □ 



From Remark [33] we have that 

Q = 1 - 2gp(0)Y(pol 



As we stated in the proof of Lemma 3.4 the right-hand side above vanishes for A = Aq, and approaches 1 
monotonically from below as A — » oo. Therefore, for some A c , we have 

4 ^ o; A) < f mx I II i r for a11 A > **■ (3 - 14) 

gp(0) + g\\Pp\\L"\Po\ 

As 2gp(Q)Y{pQ\ Aq) = 1, and Q obtains its unique minimum at Aq, we must have, in particular, that A c 
occurs to the right of Aq. The significance of this quantity will be made clear later. Briefly, the point is that, 
when A* > A c , we will show that the null space of the linearized problem is one-dimensional, whereas in 
general it may be two-dimensional. In fact, when p is taken to be a constant, one can instead take A c — Aq 
(cf. 11321 ). in which case it satisfies the implicit relation 

r° 1 

H p (p;Ac?dp = ——. (3.15) 
Jpo SP(0) 

Since H p {-\ A) is strictly monotonic as a function of A, ( |3.15| l serves as an equivalent definition of ^ c when 
p is constant. 
Next, put 

cr c := (gp(0)) 2 f (HpipUcT 1 +gP P ip)){[ P H p (s;A c ) 3 ds) dp. (3.16) 

J Po Wpo / 

The importance of this quantity is found in the following lemma. 
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Lemma 3.9 Suppose that po, cr, p, and fj are given satisfying ( |L-B[l (or size condition \\.\5\ ). Then, for a 
unique A„, there exists a solution of the form m = M(p) cos q to ( |3.3[ ). Moreover, if we assume that 

<t>o- c , (3.17) 

then A c < A*. (Note that this is an explicit statement as neither H nor p depend on cr.) 

Proof Writing ( |L-B[ ) in terms of p we find that po, p and p satisfy the Local Bifurcation Condition of [36|. 
Thus applying Lemma 3.4 through Lemma 3.9 of that paper gives the existence of a solution M of ( |3.10| l for 
a unique A*. As we have commented before, this is equivalent to solving ( |3.5| ), so m(q, p) := M(p) cos q is 
indeed the solution to ( |3.3| l we seek. 

Next, we the minimization problem ( |3.1 1[ ). For each A, the function M that attains the minimum, p(A), 
will satisfy the boundary conditions of the linearized problem ( |3.10| ) as well as the ODE: 

[HfM p } p = -^A)(H-' + gp p )M. 

By construction, therefore, p(A,) — -1. Fix A = A c , and let <p e S be defined as 

(p(p) = I a(s)~ 3 ds, po < p < 0. 

J Pa 

Note that <fi is therefore the solution of the following boundary value problem: 

(a 3 cf, p ) p =0, f OI p Q <p<0, 0(po) = O, A 3/2 cf> p (0) = gp(0)4>(0)- 
Evaluating "R(cf>; A) yields 

-gp(0)m 2 + £ a(p?<f> p (p) 2 dp 

*c) = 5 

I Po ^ 2 { a{ P"> + SPp(P)) dp 

-cr<f>(0) 2 

f po <I>(P) 2 (a(p; X) + gPp(pi) dp 
= -o-crf. 

For cr > cr c , the above identity implies fi(A c ) < - 1 . By the monotonicity of p, we conclude that A c lies to 
the left of A„ as desired. □ 

We end this subsection with a lemma characterizing the set of cr for which zero is an eigenvalue of the 
linearized problem. 

Lemma 3.10 (Zero Eigenvalue) Fix pa, p and j3. There exists at most one value of cr for which both ( |L-B| l 
holds and there exists a nontrivial (M, n) solving \3.5\ with n = and A = A*. In particular, £3 is either a 
singleton or the empty set. 

Proof Let o~ Q be given as in the statement of the lemma. Then by ( |3.5| l with n — 0, there exists M solving 

(a 3 M p ) p - g Pp (M - M(0)) = 0, M(p ) = 0, A 3/2 M p (0) = gp(0)M(0), 

for A = A*. Note that all the effects of surface tension have been removed from the equation, as cr occurs 
only as a coefficient of n in ( |3.5| l. Indeed, this problem is identical to that studied in [36], where it was 
shown that a nontrivial solution exists if and only if A t — Aq. 
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In order to prove the lemma, therefore, we must show that there is at most one value of cr for which A* 
can coincide with Aq. To make the dependence of A t on cr explicit, define 



-gp(0)<A(0) 2 - am 2 + £ a 3 <f>ldp 
IpJ 2 { a + 8P P ) 



n&A,cT) := —dp, 



for (p e S. Then, for each cr for which ( |L-B| > is satisfied, we may let A t = A t (cr) to be the smallest value of 
A > -2B mm + eo such that 

M<R((f>;A t (cr),o-) = -1. 

Clearly, /I, is continuous in cr on its domain. 

Now let <j\, a"2 be given with en < o~ 2 , ( |L-B[ l hold for crj (hence also for cr 2 .) By the definition of A*(-), 
there exists (pi, <p 2 £ »S such that 

^ij/UCcrj),^) = -l =<R(<p 2 ;Mo-2),o- 2 ). 

But, recalling the meaning of %(■; ■, •), we have 

9?(0 2 ;^(tr 2 ),Cr 2 ) < W i; ^(cr 2 ),cr 2 ) 

= Kft;i(o- 2 ),(7i) + — q- 



JL ^ (a + gp p )dp 



Jpo 

<K((Pi;A*((T 2 ),(Ti). 

The strict inequality derives from the fact that (pi(0) + 0, in light of the equation ( |3.5| l satisfied by (pi. We 
have therefore shown ^((^/^(o^Xcri) > -1. By the strict monotonicity of p, we conclude that A„{<T2) 
lies to the right of A t (<T\). It follows that A„ is a strictly increasing function of cr. Hence there is at most 
one value of cr for which A,(cr) = Aq. □ 



3.3 Proof of local bifurcation from simple eigenvalues 

In short, the previous section assures us that, if ( |L-B[ ) and size condition ( |3.17| l each hold, then, at some 
A t > -2B mm + €q, there exists a solution of the linearized problem \33\ . In order to infer the existence of 
curve of non-laminar flows local to this point, we shall make use of the classical theorem of Crandall and 
Rabinowitz. In this section and onward, we shall use N{£) to denote the null space of a linear operator £, 
and to denote its range. 

Theorem 3.11 (Crandall-Rabinowitz, [8|) Let X and Y be Banach spaces, I c R an open interval with 
A, e I. Suppose that T '. I X X ^> Y is a continuous map with the following properties: 

(i) T(A,0) = 0, for all A el; 

(ii) D\T, D2T and DiD 2 T exist and are continuous, where D, denotes the Frechet derivative with respect 
to the i-th coordinate; 

(iii) D 2 !F(/l»,0) is a Fredholm operator of index 0, in particular, the null space is one-dimensional and 
spanned by some element w*; 

(iv) DiD 2 f(A t ,0)w t t t R(D 2 T(A t ,Q)). 
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Then there exists a continuous local bifurcation curve {(A(s), w(s)) e RxX : |s| < e) with e > sufficiently 
small such that (/l(0), w(0)) = (A t ,w„), and 

{(A, w) e U : w * 0, T(A, w) = 0} = {(A(s), w(s)) e R x X : |j| < e} 

/or some neighborhood 11 of (A*, 0) in 1x1 Moreover, we have 

w(s) = w* + o(s), i« X, |j] < e. 

IfD^F exists and is continuous, then the curve is of class C 1 . 

All that is left for us now is to verify the hypotheses of the preceding theorem. As in the previous 
sections let the transformed fluid domain be 

R := [(q,p) : < q < 2n, po < p < 0}, 

with boundaries 

T := {(q,p) eR:p = 0}, B := {(q,p) e R : p - p ], 

and define 



X := (ft e C^(J!) : h = on S}, Y = Y, x 7 2 := C^ tt (7?) x C^(T). 
Let h(q, p) :- H(p) + w(q, p). Then by the full height equation, w must satisfy the following PDE 



(1 + w 2 q ){H pp + w P p) + w qq (H p + Wp) 2 - 2w q Wp q (H p + w p ) 



(3.18) 



-g(H + w- d(H) - d(w))(H p + Wpfpp + (Hp + WpY/Si-p) = in R, 

I +w 2 + (Hp + w p ) 2 (2ctk[w] + 2gp(H + w) - Q) - on T, (3.19) 

together with periodicity in q and vanishing on B. 

Conforming to the framework of [8|, we introduce a nonlinear operator 

r = (Tu Tz) : (-2Bmin + eo, co) x X -> Y 

defined, for w e X, A > -2B m ; n + eo, by 



Ti(A, w) := (1 + w 2 )(H pp + w pp ) + w ?? (H p + w p ) 2 - 2w q w pq (H p + w p ) 

-g(H + w- d(H) - d(w))(H p + w p f Pp + (Hp + w p ) 3 B(-p) 



(3.20) 



T 2 (A, w) ;= 1 + w\ + (Hp + w p ) 2 (2o-K[w] + 2gp(H + w) - Q). (3.21) 

Note that by definition of the laminar solution H(-;A), we have T(A, 0) = 0. For future reference, we 
evaluate the Frechet derivatives T\ w , Tzw at w — 0: 

r lw = d 2 p + H 2 p d 2 q + 3H 2 p B(-p)d p - 3g(H - d(H))H 2 pP pdp - gH 3 pPp (l - d) (3.22) 

r 2w = (2gpH 2 p + 2Hp(-2crd 2 + 2gpH - Q)d p ) 



(zgpJT 1 - 2A l/2 d p - 2crA~ 1 d 2 )j 



(3.23) 



In order to show that the eigenvalue at A t is simple we must characterize the null space and range of 
the linearized operator J\ v (A t ,0). This is accomplished in the following three lemmas 
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Theorem 3.12 (Null Space) For A* + Aq, the null space of f w (A*,0) is at most two-dimensional. In 
particular, it is one-dimensional provided A„ > A c . Zero is an eigenvalue if and only if A, - A$. 

Proof In the previous section we showed that M(p) cos q e N(Tl v (A t ,0)), hence the null space is at least 
one-dimensional. Using a generalization of the arguments in [321, we now demonstrate that, when A* + 
Aq, then AfC7v(/L, 0)) can be at most two-dimensional. Furthermore, when A* > A c it is at most one- 
dimensional. 

Let m e NCF w (A*,0)) be given. Then, by the evenness of in in q, we may decompose it via a cosine 
series: 



m(q,p) = ^ m„(p) cos (nq). 



n=0 



By the definition of A*, we have that m\ is in the span of M, the minimizer of the Rayleigh quotient H. 
Because m is in the null space of the linearized operator, each m n . must satisfy (|3.4|) for n — n t . 



The case when mo £ was dealt with in Lemma 3.10 Recall that it was argued that this occurs precisely 
when A t — Ao. 

Let us now consider the case A„ + Aq. Suppose that n\ , «2 an d «3 are distinct positive integers such that 
1 e {«i,«2> w 3} an d m ni ,m„ 2 ,m„ 3 £ 0. We shall prove that (i) for A > A c , each m n . is in the span of m\, 
and (ii) at most two of {m ni ,m n2 ,m n] } are linearly independent. For simplicity let us denote m n . :— 0,- for 
i = 1,2,3. 

As «; > 1, from \3.5\ we conclude that the functions </>, solve the following ODE: 



-(a 3 #)' + gp'(pt = -nfa4>, p <p<0, 

along with the boundary conditions 

0,(p o ) = 0, A^iO) = (n]o- + gp(O))0(O), i= 1,2,3- 
Multiplying the first equation by cf>j and integrating, we obtain the following identity: for any p < p < 0, 

p pp r-p 
(a 3 ^) = a 3 (f>' i cf>' j dr+ (nfa + gp')(p i (f>jdr>0. (3.24) 

J pa J pa 

When |? = 0we find from the boundary conditions that 

- n? ( f aptfjdp - 0-0,(0)0/0)) = -gp(O)0,(O)0 7 (O) + f c?^)dp + f gp'^jdp- (3.25) 
In particular, when i - j we get the useful identity: 

- n] I f acfidp - o-0,(O) 2 ) = -£ P (O)0,(O) 2 + f a\^dp + f gp'0 2 </p. (3.26) 

\Jpo I J po J Po 

Note that the quantity in parenthesis on the left-hand side need not have a sign. In some sense, this is what 
is responsible for the appearance of double bifurcation. 

Importing some notation from l32l . let us define the space H := L 2 ([po,0]) X C and endow it with 
additional structure by prescribing the indefinite inner product 

f° _ - 
[ui,U2]m'-= I au\U2dp - crb\b2, foreach = (m,,o,) e H, = 1,2. 

•Jpo 
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Note that [•, -]h consists of an infinite-dimensional positive definite term and a one-dimensional negative 
definite term. As a trivial consequence, we have that any maximal negative definite subspace of H (with 
respect to [•, -]h) is one-dimensional. 

We also observe that, taking ( |3.25| l and exchanging the roles of </>,■ and <pj, reveals 

-n 2 |J afrfadp - 0-0,(0)^(0)1 = -n) ajmpjdp - tr^(0)^(0)J . 

Recalling the definition of [•, -]h> this last identity implies immediately that, 

$M = 0, for i,y = 1,2,3,1 ^ j, 

where 0, := (</>,■, 0,(0))- In other words, the 0, are mutually orthogonal with respect to this indefinite inner 
product. 

Let Dclbe the dense set D : = {(a, b) e H : u e H 2 ([p Q , 0]), u(po) = 0, b = m(0)}. In particular, we 
have by the boundary conditions imposed in ( |3.4| l that (j), is an element of D. 
Now let A > A c be given. Then, for any u e D \ {0}, we estimate 



W = \u(0)\ z = 



I u'(p)dp 

•Jpa 

^ ri + G(p-A) r* ajpa? 2 
r° 

<4*%Po) a{p;A?u'(p?dp. 

J Pa 



The last inequality come from the estimate 1 + G > 1/2 and the identity ( |3.2[ ) evaluated at p = po. Since 
A„ > A c , it follows from above that 



\b\ 2 < , n . I , , , f a(p;i)V(p)V (3.27) 



Using <\3.21\ with 5 = </>,-, we find 

-0 



f aXftfdp - sp(O)0,(O) z - gllp'MPol > 0. 

J Pa 



From ( |3.24| > with i = j, however, we conclude that </>,(-) 2 is a strictly increasing function attaining its 
maximum value at p = 0. It follows that 



r-0 

< -sp(O)0,-(O) 2 + a\<l>' i ) 2 dp+ gp'cf> 2 dp = -n 2 [4, i ,0 l 

J Pn J Pn 



Here the last equality is simply a restatement of ( |3.26[ ). But now we have shown that each nontrivial </>, 
is negative definite with respect to [■, -]e- Since they are all mutually orthogonal, this is impossible unless 
they lie within the same one-dimensional subspace of H. Thus, for A > A c , N(f w (A, 0)) is given by the 
linear span of m\ . 
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For A < A c , the quantity appearing on the right-hand side of ( |3.26| l may be negative and so we must 
approach the analysis with greater care. Let us consider the indefinite inner product 

r° _ f° , _ 

[u\,U2\d '•= I a u x u' 2 dp + I gp UiU 2 dp - gp(0)u\(0)u 2 (0), for every u\, w 2 e D. 

J Pa Jpn 



We claim that every subspace of D that is maximal and negative definite with respect to [■,■]£> must be 
one-dimensional. Since there is an infinite dimensional negative semi-definite term in the definition of 
[•> 'Id, proving this claim is not automatic. 
Let u e D be given. Then we estimate 

r° r° I P' \ 2 

I gp'(pMpfdp= I 8P'(P)\\ u'(r)dr\ dp 

r-0 r-p 

^ I gp'(p)ip-Po) I u'irfdrdp 

J pa J p a 

= I W(rf I gp'(p)(p - p )dpdr 

Jpo Jr 

i r° 

> -xgWp'h-pl u'ipfdp. 

Jpo 



Recall, however, that by the definition of e , 

a 3 > ef > 2g\\p'\\ L ~p 2 eM > 2g\\p'\\ L ~p z . 
With this inequality we are more-or-less done, for combining it with the last inequality we get 



1 f 

[», u\ D >- a 3 (u') 2 dp - gp(0)u(0f. 

^ J Pa 



The right-hand side above clearly defines an indefinite inner product on D for which each maximal negative 
semi-definite subspace is of dimension one. The same must therefore be true of (D, [>, -] D ). This proves the 
claim. 

Now, using the notation we have developed, ( |3.25| l can be rewritten: 

- = 4>j] D , U j = 1, 2, 3. (3.28) 

Since the 0, are orthogonal with respect to [•, -]h, ( |3.28[ l says they must also be orthogonal with respect to 
[•, -]d- Appealing to the claim, we conclude that at most one of them is negative semi-definite with respect 
to [•, -]d- Without loss of generality, we may therefore assume that [c/>2, <I>2]d, [03, 4>3\d > 0. But then ( |3.28| l 
implies 02 an d 03 are each negative definite with respect to [•, -]h, which is impossible unless they are in 
the same one-dimensional subspace. For general A ? Aq, therefore, we have proved N(9\ V (A, 0)) is at most 
two-dimensional. □ 

Remark 3.13 From the arguments above it is clear that, when A* = Aq, the dimension of N{T w (A t ,0)) is 
at most three. This is simply a consequence of the fact that there can be at most two nonzero eigenvalues 
— one whose eigenfunction is negative semi-definite with respect to [•, -]h and one whose eigenfunction 
is positive definite. When p is a constant, one can leverage this observation to prove that the null space 
will be at most two-dimensional. Unfortunately, for general p, <po, the eigenfunction for n — 0, need not be 
orthogonal to the other eigenfunctions (with respect to [•, -]h-) 1=1 
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The previous lemma motivates the choice of notation in p.6| >-( |378| : E\ consists of those values of <x 
for which simple bifurcation occurs; E 2 consists of those for which we have double bifurcation; and Ej 
(if it is nonempty) contains the unique value of cr for which zero is an eigenvalue and we cannot rule out 
triple bifurcation. The next corollary makes this statement more precise. 

Corollary 3.14 Fix p ,p and (5. Denote 

E :- {cr e R + : (po,cr,p,fi) collectively satisfy ( |L-B| >}. 

Then 

(cr c , 00] c E\ and E\ U E 2 U E$ — E. 



Proof Theorem 3.12 implies that, if A* + Aq, then NCFn(A t ,0)) is either one-dimensional or two-dimensional. 
First consider the one-dimensional case, i.e. assume there exists cp e X generating the null space. Since (p 
is even in q, we may write it as a cosine series: 



(f>(q,p) = V M n (p) cos (nq). 



n=0 

Observe that, because 7\,.(/l„,O)0 = 0, we have additionally 

r w (A*, 0) (M n (p) cos (nq)) = 0, for all n > 0. 

The local bifurcation condition guarantees that M\ £ 0. But, since N(^ W (A*, 0)) is one-dimensional, it 
follows that M n = 0, for n + 1 . Hence there is a unique solution (n\ , M\ ) to ( |3.5[ ) with A = A*, n\ — 1 . We 
infer that, if cr > cr c , then cr e E\ . 



In particular, note that by Lemma 3.9 when cr > cr c , then the corresponding value of A, will lie to the 



right of A c and thus the null space is one-dimensional. In other words, (cr c , 00] cE\. 

Suppose now that A, + A§ and the null space of ^.(/l^O) is two-dimensional with generators <p\ and 
cf>2- Then, decomposing each of these into cosine series we find that 

c/>i(q,p) = Mi(p) cos (nig), <pi{q,p) = M 2 (p) cos (n 2 q), 

where (ni,M\) and («2,A^2) are distinct solutions to (|3.5|l with A — A*, n\ = 1 and n 2 + 1. In fact, by 



Lemma 3.10 we actually have n 2 > 2. By definition this means cr e E 2 . 



Finally, if A t = Aq, then in view of Lemma 3.10 we have that cr e Ej, — indeed, it must be the case that 
•£3 - {c}- Since these possibilities are exhaustive, we conclude E — E\ U E 2 U E^. □ 



Lemma 3.15 (Local Fredholm Map) For each A, the map T W (A, 0) : X — » Y is Fredholm with index 0. 

Proof This follows trivially from combining Lemma 3.5 of |I32] and Lemma 4.2 of ll36ll . □ 

Lemma 3.16 (Range) The pair (^l,S) belongs to the range of the linearized operator ^.(/l^O) if and 
only if it satisfies the orthogonality condition: 

1 ' „2,* 



J J Jla 3 fdqdp 'Sa dq = 0. 



where cf>* generates the null space ofT w (A t ,0). 
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Proof Suppose first that (Jl, S) e Y is in the range of T w (A t , 0). Then there exists some v e X such that 
J{ = :Fi w 0i*, 0)v, and £ = r 2w (A», 0)v. It follows that 

jj 3la 3 cf>*dqdq = JJ ((a 3 v p ) p + av qq - a 3 g Pp v)cf>*dqdp 

= jj ((a 3 4> P )p + af qq - a 3 gp p (f>*)vdqdp + J (a 3 v p <f>* - a 3 vf p )dq 
= J (a 3 v p (p* - a 3 vf p )dq, 

as on B, both v and 0* vanish. Now, on T we have a 3 — A 3 ^ 2 . Moreover, from the definition of cf>* we have 
on T that gpcf>* - cr0* = A 3/2 <f> p . Thus, 

i J" &a z fdq = J" (gp/T'v - / l 1/2 v i , - (rF l v qq )ty*dq 

= f ((gf»f>* -o-f qq )v-A 3 l 2 fv p )dq 

= £(A 3 %v - A 3 ' 2 fv p )dq. 

Combining this expression with the last gives necessity of the orthogonality condition. 

Sufficiency follows in the same way as Lemma 3.6 of [32|. The main point of that argument is that, 
once we have established TwiA, 0) is Fredholm of index 0, then dimension counting and necessity imply 
sufficiency. For brevity we omit the details. □ 

Finally, we must ensure that the so-called transversality or crossing condition holds. 

Lemma 3.17 (Transversality Condition) Suppose that <j>* generates the null space ofF(A t ,0). Then 'F w a.( / ^*< 0)0* ( 
7?CFwW„0)). 

Proof First we calculate the mixed Frecfhet derivatives of f at (A t , 0): 

TuJA t ,0) = -(1 + G)a- 4 d 2 q - 3(1 + G)a A B{-p)d p - 3gYa 2 p p d p 
+ 3gY(l + G)a A p p d p + 3 -g(\ + G)a- 5 Pp (l - d) 



T 2 M,0) = ( - 2g P A' 2 - A- l/2 d p + 2crA- 2 d^ 

By the previous lemma, it suffices to show that the pair CFu w (/l*, 0)0*, T2Aw(^*, 0)0*) does not satisfy the 
orthogonality condition. Equivalently, if we put 

3 := H( (i + 6)a ~ i(<p * )2 ~ 3(1 + & a ~ i p(-pw*<pp 

- 3gYa Pp ff p + 3gY(l + G)a~ l Pp ftf p + + G)a" 2 P p(0*) 2 )^ 
+ J ( - gpa~ 2 (<P*f - ^aff p + aar 2 <p* qq 4>*yq, 



our lemma will be proven if we can show 3 £ (again, because d(c/>*) = 0). We will demonstrate that, in 
fact, 3 < 0. To keep our notation concise, let 3 = 3\ + . . . + <5g, where 3j denotes the i-th term in the sum 
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above. Note that the only additional term from the 3 considered in Lemma 3.12 of [36] is Sg. By simply 
integrating by parts once, and noting the boundary terms vanish by periodicity, we find 



-crA. 



f T <t> 2 q dq<0. 



The lemma follows. 

Combining these lemmas, we now prove the main theorem for simple eigenvalues. 



Proof of Theorem 3.1 That f satisfies conditions (i)-(iii) of Theorem |3.1 1| is elementary. Likewise, 
< F ww (A*,0) exists and is continuous. Finally, conditions (iii)-(iv) follow from Lemma 3.15 and Lemma 3.17" 
respectively. 

We are therefore able to conclude the existence of a C 1 -curve C[ of non-laminar solutions to the 
height equation bifurcating from T at (H(-;A*), Q(A*)). In a neighborhood H of the bifurcation point, 
moreover, this curve and T together comprise the complete solution set and each h e C\ fl K can be 
written 

h(q,p) = H(p; A*) + eM(p) cos q + o(e), inZ. (3.29) 

Since H p (-; A) > 0, p. 29 1 implies that, possibly by shrinking U, we may assume that h p > on C[ . 
Then, in light of Lemma 2.1 there is a corresponding C 1 -curve C\ oc of solutions to the original problem, 
< | 1 . 1 [ ) — ( fTT8| > - Since C( n T — {(/l*,0)}, we are guaranteed that C\ oc contains only one laminar solution: 
that corresponding to H(-;A,). Finally, part (iii) of the Theorem statement follows from ( |3.29| l, possibly 
by further restricting to a subneighborhood of (H t , Q t ) in 14. □ 



4 Local bifurcation from double eigenvalues 

If we do not assume that cr e E\, then we cannot rule out the possibility that null space of ^(/L,0) is 
two-dimensional (or, indeed, that zero is an eigenvalue.) Were this to be the case, the Crandall-Rabinowitz 
arguments of the previous section break down and we are forced to revert to a more general methodology. 
To successfully prosecute this program requires that we impose some additional nondegeneracy condi- 
tions. Suppose that cr e £ 2 U Z3. We may therefore let + 1 be given so that there exists a nontrivial 
solution to (j33J with A = A*, n = n-i. We require 



and 



« 2 > 3 (4.1) 

* 01U102322 * ^112202211, 01111^22^02211^11 #2222 f 11 * &1 122^22, (4.2) 

where &ujj, 'Pu are given by ( |4.25| l-( |4.26| l and ( |4.14| l, respectively, for i,j= 1 , 2. In particular, we note that 
( |4.1[ ) implies cr <£ 2T3. The purpose of dictating 112 + is merely to avoid considering the zero eigenvalue 
case, which require some additional work. In fact, Wahlen was able to give a fairly complete analysis of 
the local bifurcation when this does occur, though the argument is quite subtle (cf. II32I ). The reason for 
avoiding n 2 — 2 is less obvious, but results from certain symmetries in the underlying problem that shall 
become clear later (in particular, see Lemma [43] ) 

Under these assumptions, the eventual product of our efforts in this section will be the proof of The- 
orem 1.4 We shall accomplish this in several stages, according to the Lyapunov-Schmidt procedure. The 
techniques here are not new, though, as we shall see, there is a certain degree of degeneracy in the problem 
that will require some care. 
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Step I: Reduction to finite-dimensional problem. We begin by recalling a weighted inner-product on Y 
introduced in the previous section. For all pairs ( Jl,, S,) e Y, i — 1,2, let 

muS 1 ),(M 2 ,S2)) Y := ff a 3 (-;A,)Jl l (-)Jl 2 (-)dqdp 

1 r (4.3) 
+ 2 J a\-;A t )Bi{-)B 2 {-)dq. 



Denote Xa := N(T w (A t ,0)). Then, as we have seen in Theorem 3.12 Xa will be at most dimension two. 
Suppose that it is spanned by <f>\ and (f> 2 , where <f>i(q,p) = M,(p) cos / = 1,2. Of course, the most 
interesting case given our previous analysis, is when n\ = 1. But there is little to gain at this stage in 
making such a restriction. 

It follows directly from Lemma 3.16 that the range of Tw(A t , 0) is given by the orthogonal complement 
(with respect to (•, -)y) of the linear span of [<f>i,(f>2}- So, if we now define 

X B := {heX: ((h,h\ T ),(<f>u4>i\T)) Y = 0, i = 1,2} , 

Y A := R^l,S) e F : &\ T = fl, e spanKfc.fclrMfe.fclr)}}, 

Y B := {{ft,®) eY:3l\ T =&, ((&,&), (fc, fc|r))y = 0, * = 1,2}, 

then X = Xa © X B , Y - Y A ® Y B . Finally, let A x : X — » Xa, Sx : X — > X B denote the projections onto X A 
and X B , respectively, and similarly for Ay : Y — » Fa, 5y : Y — > 7^. 

As always in bifurcation theory, we are trying to find all solutions to the problem 

lF(A,w) = 0, (4.4) 

in a neighborhood of (A„,0) e R X X. Applying the projections above, we see that this is equivalent to 
solving 

A Y T(A, A x w + B x w) = 0, (4.5) 
B r T(A, A x w + B x w) = 0, (4.6) 

locally near (/t„,0). 

Consider ( |4.6| l . Merely by unravelling the definitions above, we have A Y T K (A* , 0)B X = and therefore, 
for any w eX, 

B Y T w (A*,0)w - r w (A*,0)B x w = B Y r w (A,,0)(l - B x )w = B Y T w (A t ,0)A x w = 0. 
This proves that Tw(A,,0) commutes with projections onto the infinite dimension spaces in the sense that 

r w (A,,0)B x = B r r w (A,0). (4.7) 

Now we note that, by construction, T w (A t ,0)B x is an isomorphism from X onto Y B - 'R(T K {A t , 0)). 
Thus, in view of \A.1\ , the Implicit Function Theorem allows us to solve for B x w as a function of A x w 
and A near (A*,0). In particular, if we write A x w = gi<p\ + £ 2 <p2, for some £1,^2 £ R, then there exists 
a neighborhood U c R 2 X R of (0, 0, A,) and a smooth £ : t/ — » 7 B such that, in a sufficiently small 
neighorhood of (A t ,0) inRxX, all solutions of ( |4.6| are given by {(A, %\<pi +fe02+£(£i,fe; X)) : (^\,^ 2 ,A) e 
£/}• Moreover, it is easy to see that £(■, •; A) = 0(|^| 2 ), in this neighborhood. 

Taking these solutions and inserting them into ( |4.5| l, we infer that ( |4.4[ l reduces to the following finite- 
dimensional problem: Find all 6 R 2 , 1 e R in a sufficiently small neighborhood C/ of (0,0, /l.) 
such that 

8(fi,&;A) = 0, (4.8) 
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where the bifurcation function 93 : C/ — » R 2 is defined by 

S8i(f 1 ,&U):=(0 < ,mfi0i+&02+^i,&U))) y , i = 1,2. (4.9) 

To simplify this expression, we repeatedly Taylor expand T . In doing so, we shall let 

(f w aO)^^>eX(I,}0 

denote the second Frechet derivative of T in the <p { then <f>j direction, j = 1,2. With this notation, 
expanding \A.9\ around (/L, 0) yields, for j, fc, ^ = 1,2, 

93 ; = (A - /UJfyJFy + tj{ k ijk + tMc&iju + 0(\A - A,\ z \® + \A- A t \\® 4 + |£| 5 ), (4.10) 

where we are using summation convention and 

Vij ■= \ (<Pi, TM,^i) Y , i, 7=1,2, (4.1 1) 

*ijk := ~ {<Pu (r ww U*,0)<pj, <p k )) Y , i, j,k = 1,2. (4.12) 

0yw:=g(0i. ((r www (A„0)4> j ,(p k ),<f> { )) Y , i,j,k,£= 1,2. (4.13) 



iS'fe/) 77: Computation of coefficients. Clearly, before we can begin to characterize the solutions of ( |4.8| l, 
we must first determine (P,-^ and 6> (JJ tf more precisely. As one might expect, this will be an elementary 
but involved computation. What will simplify things immensely, however, is the special structure of <p\ 
and <p2- That is, we know c/>i(q,p) = Mi(p) cos (ntq), where M, is the solution of \3.5\ with n = for 
i — 1,2. This will aid us in evaluating the integrals in (•, )y. 

Also, observe that multiplying ([33} by M, and integrating by parts one obtains, for p € (po, 0], i — 1,2, 

rp rp 
MiM'i = a 3 a 3 (M-) 2 dr + a 3 {an} + gp p )Atfdr. 

Jpo Jpo 

As the right-hand side above is strictly positive, we surmise that M, and M ( ' do not vanish on (po, 0) and 
are of the same sign. By linearity of the equation, we may simply take them both to be positive. 

Remark 4.1 In general, for a Sturm-Liouville problem one has that the only eigenfunction without nodes 
is the ground state. We have argued, however, that M,, M'. are both nonvanishing for i — 1,2. This seeming 
contradiction is precisely due to the presence of the surface tension term on the boundary, which alters the 
orthogonality relation between the eigenfunctions. □ 

Lemma 4.2 (Linear Terms) Tlte matrix fjj defined by is diagonal and negative definite. 

Proof We shall start by explicitly computing ¥u. Using the expressions for T\ W (A*, 0) found in the previ- 
ous section, we see that, for i, j = 1,2, 



f n)c 

•J Dn 



■f 

'J Po 



a-\l+G)M i M j dp-3 {\ +G)a 1 pM i M' j dp-?>g Yap p MM'jdp 



Jpo 



+ 3g [ Y(l+G)a- l p p M i M' j dp+ : -g j (1 + G)a i p p M i M j 

Jpo ^ Jpo 

1 ^ r 2n 

- - (2gpA~ l MjMj + A V2 MjM'j + 2o-A~ l n 2 ^ J cos (mq) cos (njq)dq. 
Thus, fij = for i + j. Arguing as in Lemma 3.17 moreover, we can show that f„ < 0, = 1,2. 



(4.14) 
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Let us now attempt to characterize <Pijk- Unfortunately, the conclusions we will be able to draw will be 
far less satisfying than for the previous result. 

Lemma 4.3 (Quadratic Terms) Let <L>ij k be defined as in \A.\2\ then, if nj + 2n\, we have @i jk = 0, for 
i, j,k — 1,2. On the other hand, if n-i -2n\, then necessarily, <P\\\, 0222, ^221 > ^212* ^122 = 0. In general, 
0\\2 and <P\22 ( = &211) are given by equations {X.\9\ and ( |4.18| l, respectively. 

Proof Recycling notation slightly, for the purpose of the next computation we let and be given ele- 
ments of X with d(4>) — d(i/r) = 0. This restriction is natural as we shall eventually be applying this to 
01, 02, each of which have this property. Then, 

CFnv W (A, w)<p, 0) = 2w q i// q (p pp + 2(H PP + w P p)<f>qi!f q + 2w q i/t pp <t> q 

+2(H p + W p )\fi p <p qq + 2{H p + Wp)\jf qq <pp + 2w qq <j/ p <p p 

-2w q ifr pq cp p - 2w qp \l/ q <t> p - 2(H P + Wp)\jj q ^ pq - 2w q i]jp(pp q 
-2w pq i/f p (f> q - 2(H P + Wp)i// pq (/) q - 3gp p (H p + w p fft p (/> 
-3gp p (H p + Wpftfnfip + 6(H P + w p )itf p (p p p(-p) 
-6gp p (H p + w p )(H + w - d(H) - d(w))ilf p <pp, 



and 



where 



{TiwwiA, w)cf>, 0) = 2cf> q i// q + 2(p p if/ p (2ctk[w] + 2gp(H + w)-Q) 
+2(H p + w p )<fr p (W[w]0 + 2gp4>) 
+2(H P + w p )(p p (2ctk'[w]iI/ + 2gpiff) 
+2o-(H p +w p ) 2 (K"[w]<f>,4,), 

, 3(c/> qq iff q + ^ qq ) \5w q w qq 

(k M0,0) = — 2 5/2 zwi ^v 

(1 + w q y /z (1 + w q )'' z 



Evaluating each of these at (/!*, 0), we arrive at the following 

(Ti ww (A t , 0)0, 0) = 2H pp (f> q iff q + 2H p i[f p (f> qq + 2H p ^i qq (p p - 2H p i// q <p pq 

- 2Hp4f qp <t> q - 3g Pp H 2 ^ p ^ - 3g Pp H 2 p W p (4. 15) 

- 6gp p H p (H - d(HM p cf> p + 6Hp^ p p(-p) 
(f2 W w(^, 0)0, 0) = 20,0, + 2(f> p iff p (2gpH -Q) + 2H p il/ p (2gp(/> - 2cr0 M ) 



+ 2H p <f) p (2gpif/ - 2cnj/ qq ) + do-H 2 {<p qq ^ q + (/> q i// qq ). 



(4.16) 



iprpyt-SFV ^Vqq) > i-L p yy qq r q ' y q v qq ; 

Note that, as the above expressions make clear, the regularity of T WK guarantees that {T mv (A, 0)0, 0) = 
CF WW (A., 0)0, 0). So, recalling the definition of 0,j k in ( |4.12[ ), we see immediately that 0,j k = 0/ k j, for 
i,j,k = 1,2. 

To keep things manageable we will evaluate the inner-product in ( |4.12[ ) in two steps. First we compute 
weighted the L 2 inner-product over R: 



C (2 > 

ijk 



(« 3 0,-,(^ H ,,(i„O)0 7 ,0 it » L2 = eg [n j n k (a 3 H pp M i ,MjM k ) Ll -2n j n k (a 2 M i ,(MjM k y) L2 \ 

2nj (a 2 M t , MjM' k ) L2 + 2n\ (a 2 M /; M)M k ) Ll 
- 3g (ap p Mi, (MjM k )') L2 - 6 (a 2 pM h M'jM,)^ 
+ 6g (a 2 p p (H - d{H))M h M'jM' k ) L2 
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where, for i, j,k = 1,2, we put 

cos (n,#) sin (rijq) sin (n k q)dq, 

X 2 n 
cos (riiq) cos (njq) cos (rikq)dq. 

Similarly, evaluating the weighted L 2 inner-product on T yields, 

\ (a 2 <f>i, (r 2ww (A t , 0)<Pj, <Pk)) L2{T) = Cg [n i n k A t M i M j M t + 3n 2 n k <rM i M j M k ] 



+ C (2) 



ijk [Al' 2 (2gp + 2crn 2 k )M i M' j M k 
+ A l J 2 (2gp + lo-rfyMiMjMl 
- A 2 MiM'jM' k 
+ 3C ( f k o-n j n 2 k M i M j M k , 



where 

- 2 n 



Xz/r 
cos (iiiq) cos (njq) sin (n k q)dq, 



Xz/r 
cos (ntq) sin (njq) cos (n k q)dq. 



On both the boundary and interior, therefore, we will have massive cancellation due to the integrals of 
the various trigonometric terms. For instance, when = j = k, all of these integrals vanish, which implies 
0m = 0222 = 0. More strikingly, unless n 2 = 2«i, every term will vanish and (J * = 0, for i,j,k =1,2. 
Suppose that we do have this relationship between n 2 and m, then, after the dust clears, we are left finally 
with the following: 

0111, 0222, 0221, 0212, 0122 = 0, (4.17) 

^0211 = -n 2 (a 3 H pp M 2 ,M 2 ) L2 -3g(ap p M 2 ,(M 2 y) L2 

- 6g (a 2 p p (H - d{H))M 2 , {M\ ) 2 ) l1 + 6 (a 2 (JM 2 , {M\ ) 2 ) l1 (4. 18) 

+ {AA\ l2 (2gp + 2crn 2 )M 2 M l M[ - A„n 2 M 2 M 2 - A 2 t M 2 (M[) 2 ^ , 



and 



^0112, ^0121 = n in2 (a 3 H pp M 2 ^ 2 ) L2 - Inl^M 2 ^!,)^ - n 2 2 (a 2 M 2 AM 2 Y) 



L- 



3g (ap p M x ,(M x M 2 )') L2 - 3g (a 2 p p (H - d(H))M' 2 , (M 2 )') l2 

+ 3 (a 2 l3M' 2 , (M 2 Y) l2 - 2n x n 2 (a 2 M x ,(M x M 2 )') ^ (4 19) 

+ (A t nm 2 M 2 M 2 + A x J 2 (2gp + 2crn 2 2 )M 1 M' 1 M 2 
+ Al /2 (2gp + 2<rn 2 )M 2 M' 2 - A 2 t M x M' x M' 2 ) . 



This completes the lemma. 
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Given the high degree of degeneracy at the quadratic level, in order to characterize the complete solu- 
tion set of the bifurcation equation we forced to consider the cubic terms. These we compute in the next 
lemma. 

Lemma 4.4 (Cubic Terms) Let be given as in ( |4.13| l. Then 

01112, 01121, ^1211, 02111, 02221, 02212, 02122, 01222 = 0. 

The remaining coefficients are generally nonvanishing; ©an is given by equation \A.25\ , while &ujj is given 
by ( |4.26| l, for i,j- 1,2 and i + j. 

Proof As in the previous lemma, we begin with an explicit calculation. Differentiating once more in w the 
expressions for T w „(k, w) found in Lemma 4.3 we find, for 0, iff, e X with d(c/>) = d{if/) = d(9) = 0, 

((TuvwwW, W)<p, iff) , 8) = 2<p pp lff q Gq + 24> q lff q 6 p p + 2<p q lff pp 6 q + 2(f> qq if/ p 9 p 

+2(ppif/ qq 6 p + 2cppiffpG qq - 2(p pq iffp6 q - 2(f> pq iff q p 
-2(f>piff pq 8 q - 2cf> q - 2cf> p iff pq 8 q - 2(p p iff q pq - 2(f> q iff p G pq 
-6gp p (Hp + w p )((f>pifrp8 + (f>piff0p + (f>iff p 0p) 
-6gp p (H + w- d(H) - d(w))(f>piffp6 p + 6j3(-p)(p p iff p p , 

and, on the boundary, 

((T 2www {A, w)cf>, iff), 6) = 2<(> p i/f p (2o-K'[wW + 2gp9) + 2iff p 9p(2o- K '[w]cp + 2gpcf>) 

+2cf>p0 p (2o-K'[w]iff + 2gpifj) + Acr{H p + w p )<f> p (K"{w\6, <A> 
+4a(H p + w p )iffp(K"[w]0, <f>) + 4cr(H p + W p)9p( K "[w]<p, if/) 
+2o-(Hp + Wp) 2 {W"[w](f>,iff),0), 



where 



{{k [w]cf>, iff), 9) = -15 +105 . 

(i + w i a )'i i (l + w q r /z 



(4.20) 



Note that this last line implies k"'[0] = 0. 

Evaluating each of these at (A*,0), we compute 

(CFiwwwU*, 0)0, iff), 6) = 2cf>ppiff q B q + 2<f> q iff q 6p P + 2(f> q iffp P q + 2<p qq iff p 6 p 

+ 2<t> p iff qq p + 2<f> p iff p Q qq - 2(f> pq iffp0 q - 2(f> pq iff q 0p 
- 2<j) p iff pq 9 q - 2cf> q - 2tp p iff pq 6 q - 2cf>piff q pq 
- 6gppH p ((ppi]fp0 + (f>piff0p + <fnff p 9 p ) 
- 6gp p (H - d{H))<f) p iff p e p 

+ 6/3(-p)(f>piff p p - 2(f) q ifjp0p q , 

({T 2www (A t ,0)cf>, iff), 6) = 2(Ppiffp(-2o-0 qq + 2gp0) + 2iff p 9 p (-2o-cf> qq + 2gp<f>) 

+ 2(f>p0p(-2o-iff qq + 2gpifj) 

+ 12crHp(f>p(i/f qq q + iff q qq ) (4.21) 
+ \2o-Hpif/p((f) qq q + <f> q G qq ) 

+ \2o-Hp0 p {(f) qq iff q + <p q iff qq ). 
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We remark again that, owing to the regularity oif — and as ( |4.20| >-( |4.21| l make explicit — we have 
Oijiit - @im(j)&(k)m(e)> where vj is any permutation on {1,2}, and i, j, k,€ = 1,2. 
To compute the inner product (■, -)y, we first calculate that, by ( |4.20[ ), 

- (a 3 (/>i, ((T lwww U*,0)<Pj, <f> k ) , <P()) L2(Kj = n k n t C^ u (a 3 MM';, M k M e ) L2 

-(a 3 MM)AM k M e )') L2 

+ n j n e C%[(a'M i M' k \Mj,M c ) Ll 

-(a 3 MM k ,(M h M?j) L , 
(a 3 MiM' { ',MjM k ) L2 
(a 3 MiM' { ,(MjM k )') 

n) (a 3 MiMj,M' k M' ( ) 
- n\ {a 3 MM), M k M' { ) 
■ n] [a 3 MM), M' k M { ) L2 
3g(p p a 2 MM),(M k M e y) L2 
3g( Pp a 2 MMj,M' k M' e ) L2 

+ 3(a 6 H pp MM),M' k M' ( ) L1 



+ njn k C (3 l ( 



r (4) 



+ / 
+ / 

+ : 
+ : 



to 



where the constants above are defined by 

-In 



ijke 



I cos cos (njq) sin (n k q) sin (ticq)dq, 
Jo 

pin 

I cos (ritq) sin (n jq) cos (n k q) sin (ri(q)dq, 
Jo 

pin 

I cos (n,fl) sin (n t q) sin (n k q) cos (ri(q)dq, 
Jo 

pin 

I cos (n,q) cos (njq) cos (n k q) cos (ti[q)dq. 
Jo 

Likewise, on the top we find 

^ (a%, ((r 2www a*, 0)0;, fa) , <f>c)) = Cg f (a 2 MM)M'M' t ) 



<(3) 

ijkt 



>(4) 
ijkt 



(4.22) 



(4.23) 



(4.24) 



Here we have used the boundary condition on T to simplify. 

As we saw in the analysis of the quadratic terms, the symmetries of the underlying equation will be 
expressed through the quantities in ( |4.23[ ). In particular, we can check that, for i, j, k, ( with an odd number 
of ones and twos, all of the coefficients vanish. Hence 0im — in this case. If there are an even number, 
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then by the permutation property observed above, it suffices to consider two types of indices: (i) 6*,,,, and 
(ii) ©njj, for i, j = 1,2 with + j. A truly elementary calculation tells us 



1 r (1) = - c 



(4) 



3n (4) _ 7T 

4 ' 2' 



while C {2) ., C <2) ., C (3) , C (3) . = 0, for i, j = 1,2 with i + j. Using the information we obtained from M.221- 
( |4.24[ ) we find that, for such j, 



and 



\e m = 1„ 2 [(« 3 m 3 ,m;') l2 -2( fl 3 Mf,(M;) 2 ) L2 ] 

- ^ [« 2 ( fl 3 M 2 , (M;) 2 ) l2 + 3g (p p fl 2 M 2 , (M,') 2 ) . 
+ (a 6 H pp Mi, (M,') 3 ) i2 j + | ( fl 2 M,- (M,') 3 ) 



i<9 l07 = \n) [( fl 3 M,M,",M 2 ) A2 - 2( C( 3 M 1 .M;,M j M;) l2 ] 

" ^K( fl3M ''(^) 2 )^ +2n j( a3M - M i> M j M 'i) L , 

+ 6g (pp^MiM'i, MjM'j) L2 + 3g (p p a 2 M 2 , («y) 2 ) 

+ 3 (c^HppMiM'i, (Mf)^ ] + \ (a 2 MM'; (m'^ 

This completes the proof. 



(4.25) 



(4.26) 



□ 



Step III: Control of higher order terms. The next step is to use the scaling in ( |4.10| l to show that all 
solutions (locally) must lie in a particular subregion of U . In the case of mixed solutions, this effort is 
complicated by the fact that 3>y# is degenerate in the £,\ -direction. Consequently, in order to apply the 
standard theory we must restrict ourselves to working in sets that lie in complements of cusps surrounding 
the £i-axis. 

Lemma 4.5 (Scaling of Solutions) (a) Suppose that 

0112,0121,0211*0. (4.27) 

For each C > 0, there is a constant K > and open subregion V c U with (0,0, /!*) e V such that, if 
93(£, A) = Ofor some (£, A) e V\{(£, A) : lfil < QfrD, then (£, A) satisfies the scaling relation |£| < K\A-A„\. 

(b) Suppose instead that 

0112,0121,0211=0, + &i m 6>2222 * ®\ 12202211- (4.28) 

Then there exists a constant K > and and open subset VofU with (0, 0, A*) e V such that, ;/23(£, A) = 
for some (£ A) e V, then |£| < K\A - /l„| 1/2 . 
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Proof First consider the situation in part (a). By contradiction suppose that, for some C > 0, there exists a 
sequence {{g n \ A„)} c V \ [(£, A) : |£| < C|£>l) such that A n \ A*, g n) as n oo, and 



1^,-^1 



as yi — * co. 



We may therefore write A n - A, - p„|£ <n) |, where p„ — > as n — » oo. Then the equation for 25i in ( |4. 10| > 
becomes, for all n > 1, 

= g> \^\ Pn W n + 2^ ) 4" ) 0„ 2 + o(|^f ) . 

Diving through by |£ (n) | 2 and taking n — » oo, each term above but the second clearly vanishes. So, in 
order to produce a contradiction we will prove that the second term is bounded uniformly away from 
zero. By construction we know |£j n) | > C|£ n) l, for n > 1. Suppose that there exists a C > such that 

C ,£«|2 < |£«)|2_ Theni 



'1 ^2 



|£(«)|2 



12 



c 



1 +c 



Ifml >0. 



Thus, we are we done if we can produce such a constant C 

Suppose no such C" exists. Then, possibly by passing to a subsequence, it must be the case that 
l^'Vl^i"'! — » as n — > oo. By the equation for ©2 we have, Vn > 1, 



= ^ |^ ( "V»^2 + i^f'P 211 + 0(|^"f) . 



As before, we may divide through by |£ (n) | 2 and take « — » 00 to find that the first term and the higher order 
terms vanish in the limit. But, 



lim 



(n)|2 



21 1 



I* 



211 



> 0, 



so we have a contradiction. This completes part (a). 

To prove (b), we shall employ a similar blow-up argument. Now, however, as the coefficients @im are 
assumed to be non-degenerate, we will not be restricted to working in the complements of cusps. 

By contradiction suppose that there exists a sequence A n )} c U with A n \ A„, t* n) —* 0, and 



as n — » 00. 



We may therefore write A„ — A, = p n |£ (,l) | 2 , where p„ — » as n — » 00. Then the bifurcation equation 
becomes, for « > 1, / = 1,2, 



<n)|4\ 



(4.29) 



Dividing through by |£ (n) | 3 and taking n^oowe see that the first term and the higher order terms vanish 
in the limit. On the other hand, by hypothesis we have, for i = 1 , 2, 



"i IT °W 



> inf If^eywl > 0. 



Hence the third term is bounded above and strictly away from zero. Having arrived at a contradiction, we 
conclude the lemma holds. □ 
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Remark 4.6 Recall that, when we prove Theorem |1.4| we will be assuming ( |4. 1 [ > — ( [472| > with n\ = 1, which 
together imply ( |4.28| l. Hence, we expect the conclusions of part (b) to hold. □ 



Fix C > and define the set W as follows 

y\{(g,^)et/:|fr|<C|ftl},iffl37 



W : 



if (4.28 



holds, 
holds, 



(4.30) 



where V is as in Lemma 4.5 Let W ± C W denote the sets {(£, A) e W : ±(A - A t ) > 0}. Then, working in 
W ± , we are justified in letting e - e(£, A) > 0, 9 — 9(£,, A) e K 2 be defined by the relations 



? = ±e0, 



\A-A t \, if (to? 



\A-A„\ 1/2 , if (|438 



holds, 
holds. 



Note that Lemma 4.5 



guarantees 



\L"(Y) 



following system of equations for (9, e) 



is bounded. Rewriting ( |4.10[ ) in terms of e and 9 we obtain the 



= 23f(0, e) 



±e j ¥ ij + e j e k & m + o{e), if 14.27 



±0;^- + Gfi k G e ijH + 0(e), if g38 



holds, 
holds, 



(4.31) 



for 1, 7, fc, ^ = 1,2. 

We wish to say that solutions to ( |4.31| l with e = are, in some sense, the only solutions to ( |4.10| l in the 
set W. To reach this conclusion requires taming the higher order terms, represented by the 0(e) in ( |4.31| l. 
We shall do this by proving the following lemma. 

Lemma 4.7 (Regular Value) Let W c U be as in (|4.30| l. Then zero is a regular value of 9 t-> SB ± (0, 0) in 
W, provided that either (a) \A.21\ holds, or (b) ( |4.28| l holds and we have, additionally, that 



01111 2*22 * 02211^11, 6>2222 ( fll * Owil^H- 



(4.32) 



Proof Assume first that hypothesis (a) holds and say that, for some 9 + in W ± , we have 23 (9,0) = 0. 



we 



Then, using our knowledge of the specific forms of 0,^ and f, ; garnered from Lemma 4.3 -Lemma 4.4 
see that, from the equation for 23*, 

±¥ n + 20 2 <2> m = O. 

Here we have used the fact that #i # 0, since 9 e W and is nonzero. But, the Jacobian of 9 t-> 23* (6, 0) is 
easily computed to be 

d& \ ±Wn it*2n n2 2 ±S 2 12 ) = ±(± ^ n + 20 ^n 2 )^22 - 4^0„ 2 cf 211 . (4.33) 

We have already shown that the first term on the right-hand side vanishes, and, by hypothesis, the second 
term is strictly nonzero. On the other hand, in the case where 9-0, ( |4.33| l shows us the Jacobian is given 
by W\ \ W22 > 0. Hence zero is a regular value in W. 

Next consider what happens when the hypotheses of (b) are satisfied. Then, using the information on 
the cubic terms obtained in Lemma |4~4l we see that 

~ ±(0 O) = /±0iFii + 0i02<9i 122 + 0j<9„n 

\ ±02^22 + ^1^26*2211 + ^2^*2222 
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Taking 8\ = we have automatically that 23* = 0. To get nontrivial solutions, therefore, requires that 82 
satisfies 

± ¥ 22 + Q\&222i = 0. (4.34) 

Note this implies that, depending on the sign of 6*22221 either 93 + (0, #2,0) = has two solutions and 
23 (0, 02, 0) =0 has none, or the reverse is true. In any event, evaluating the Jacobian of 23* at these points 
yields 

*•(*'" T ,m ^ + V J = 

Thus the Jacobian vanishes if and only if +fn + (%&\\ zi = 0. Of course, if 6*1122 = 0, this never occurs 
and we are done. Otherwise, combining this with ( |4 34[ ) leads to 

„. ^11 _.2 = _.^22 



6*1122 6*2222' 

which contradicts ( |4.32[ ). (Note that 6*1111,6*2222 ^ by the nondegeneracy condition ( |4.28| l.) Thus all 
zeros of 8 t-> 23 ± (£ i , 0) in W ± with 0\ — are simple. Due to the symmetry of 23, we can simply reverse the 
roles of one and two above to get a similar result for the case when 62 = 0, 6\ + 0. 

Finally let us consider solutions for which neither 0\ nor 62 vanish. Then simplifying we find 

/±F 11 +0?0 11 22 + 0?0im\ r4V? . 

\±W 22 + e\&22U +9*02222 J' K ' 

so that the Jacobian evaluated at such points comes to 



det 



30 



= (±W n + ^0 1122 + 3^6» lln )(±<? 22 + ^6*2211 + 3^ 6*2222) 

(9,0) 

-4^6*112202211 

= A9 X 9 2 (6*11116*2222 ~ 6*11226*2211) • 
But this quantity is nonvanishing for + by ( |4.28[ ). This completes the lemma. 



(4.36) 



Examining the proof of Lemma |4.7| we can count the number of nontrivial roots of the map h» 
23 ± (0, 0) in W under various assumptions. First suppose hypothesis (b) (i.e., ( |4.28| l and ( |4.32| l hold.) Then, 
by ( |4.34| ), taking 8\ — and <P 2222 > 0, there are no nontrivial solutions to 23~(0, # 2 ,0) = 0, and two 
nontrivial solutions to 23 + (0, 62, 0) = 0. If <P 2222 < 0, then there are no nontrivial solutions to 23 + (0, # 2 , 0) = 
and two for 23 (0, 02, 0) = 0. The same holds true when we take 02 — and search for nontrivial solutions. 
Thus there are a total of four possible nontrivial solutions where precisely one of 0\ and 82 is zero. 

On the other hand, ( |4.35| l shows us that if neither 0\ nor 02 vanish, then (0^, 8%) solve the linear system 



Note that A is invertible by assumption ( |4.28| l. Therefore, one of 23 + , 23 has no such solutions with 9\, 82 
0, while the other has four. In total, we find that there are eight nontrivial solutions of 23* (0, 0) = in W: 

(fll =0, 2 = ±| ^22 / 6*22221 1/2 ), (#1 = ±|^l/6*„ii| 1/2 , 2 =0), 

\ (4.37) 

(01 = 16*2222^11 - 6*2211 y 22 | 1/2 , 2 = ±|A|- 1 |0 lm ^ 22 - 6*1122 fll| 1/2 . 
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Suppose now that the situation in (a) occurs, that is, 0m, ^121, ^211 + 0. Then if we let B\ — 0, we 
find that 62 must be zero as well (this follows both from the bifurcation equation and, at a more basic level, 
from the definition of W.) Similarly, there are no nontrivial solutions with 62 — 0. However, if we suppose 
that 6\ , d 2 * 0, then we can solve 93f (0, e) for 2 to find 



hence, 93* (0, 0) = if and only if 



2 



0,W: 



20 



1 12 



2*72 



22 



&2U 20H20221 

Since W\ \ and ^22 are of the same sign, we see that the above equation has two solutions or zero, depending 
on whether or not 0m an d 0211 have the same sign. If 0m ^211 < 0, therefore, we find that there are no 
nontrivial solutions to the bifurcation equation, while if 0ii2 02ii > then there are two nontrivial roots 
of SB* (0,0) = 0: 



0! 



^11*72 



201120 



0, 



Pi 



112*721 



20 



112 



01 



20,120 



112*721 



20 



112 



Putting these observation together with Lemma 4.5 and Lemma 4.7 we obtain the following result 
characterizing the solution set of the bifurcation equation under various assumptions on 0,^,, &iju . 




Fig. 4.1 A typical bifurcation diagram. In this case we have <Pm, ^221 = 0, ®nn, ©2222 > and 0im02222 < 81xa.813.n- 
The two solid line pitchforks on the right consist of pure solutions, while the two dashed pitchforks on the left are mixed 
solution curves. The shaded regions are W + and W~. 



Theorem 4.8 (Bifurcation structure at double eigenvalues) Suppose that at A* > -2S„„„ + eo the null 
space off w (A f ,0) is spanned by <pj = M,(p) cos («,<?), i — 1,2, where < n\ < «2- Let <Pij, 0,jk be defined 
as above, for i, j,k — 1,2. 

(a) If <Pjj and &iju satisfy ( |4.28| l and ( |4.32| l, then there exist four C l -curves, {C[ 0G ,}-_i> of non-laminar 
solutions of the height equation bifurcating from T at A t . In a sufficiently small neighborhood of 
(Q(A t ),H f ) in X, moreover, these constitute all nontrivial solutions. Two of these curves consist of 
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pure solutions — one of minimal period 2n/n\, one of minimal period 2n/n2 — and two are of mixed 
solutions. 

(b) If instead we have that \A.21\ holds, then there exist a C 1 -curve of minimal period 2n/n2 solutions to 
the height equation bifurcating from T at A,. If additionally we have 0n2^2ii > 0, then there are 
two C 1 -curves of mixed-type, non-laminar solutions bifurcating from T at A t . In a sufficiently small 



neighborhood of the kind described in Lemma 4.5 these comprise all nontrivial solutions. 



Proof Consider the scenario in (a). As we have seen in the preceding lemmas, this means that the bifur- 



cation equation is equivalent to solving ( |4.31| ) near (9 = 0, e = 0). In light of Lemma 4.7 moreover, the 
Implicit Function theorem tells us that to each solution of 53* (9, 0) = 0, there corresponds a curve of solu- 
tions to the bifurcation equation. The arguments of the preceding paragraphs tell us that there are precisely 
four of these, two of mixed type (i.e. with 9\,9% + 0) and two of pure type — one with (9\ = 0, #2 + 0) 
and one with (9\ +0,92 - 0). The union of these curves gives the complete solution set of the bifurcation 
equation in set in a sufficiently small neighborhood of {A„, 0, 0). 

Next let us look at the situation in (b). Fix C > and let W be defined as before. Working within W we 
know that there are no pure solutions, but there are two of mixed type if and only if @m and &2i\ share 



the same sign. If this does occur, then as we argued in the previous paragraph, Lemma 4.7 and the Implicit 
Function theorem together give the existence of three curves of mixed solutions. 

However, these clearly do not give the entire solution set — Indeed, we know there must be a curve 
of purely 2w/«2-periodic solutions. This follows, as in 1321 . by noting that restricting the solution space 
X to 2n/n2 -periodic solutions reduces to one the dimension of the null space of t F w {A f , 0). Hence the 
existence of such a curve follows from a routine application of Crandall-Rabinowitz. Note that this would 
correspond to taking £1 = 0, and thus lies outside of W no matter the choice of C. □ 

We are finally ready to assemble the results of the previous steps into a proof of the main theorem. 
This amounts to nothing more than using the sufficient conditions derived in Lemma |43] and Lemma 4.4 in 



conjunction with Theorem 4.8 



Proof of Theorem 1.4 Recall that, by hypothesis, we ha ve (|4.1[ ) and \A.2\ . Since we are taking n\ — 1 



and «2 2: 3, by Lemma 4.2 



4.8 



this in turn implies ( |4.28| l and ( |4.32 l each hold. Applying part (a) of Theorem 



therefore, we conclude that there exists four C -curves of solutions, {C' lloc }'J =l , of the reformulated 



problem bifurcating from T at (Q,,H*). Likewise, by Lemma 2.1 there exists four C'-curves, {C !j i o}- = i' 
of solutions to the original problem (jTTTJ— (JTT7J . All that remains to prove are the various nodal properties 
of these solutions. 

Property (i) follows from the symmetries of the height equation, as we have argued before. Now, 
without loss of generality, we may suppose that each (Q, h) e C' jloc can be written 

h(q, p) = H(p; A) + ecf>j(p) cos + o(e), in X, 

for i = 1 , 2. In other words, we are choosing the labeling so that C\ , and C' 2 loc are the curves of pure 
solutions, while C' 3 { and C\ loc are the mixed solutions. By taking e sufficiently small (or equivalently, 
by restricting ourselves to a sufficiently small neighborhood of the bifurcation point), we have that for 
(Q, h) 6 C', , h has minimal period 2^/«, in q. Moreover, since 

h q {q, p) = -em<t> p (p) sin (mq) + o(e), in X, 

we see that, again for e > sufficiently small, h q {q, p) has the opposite sign of sin («,g) for each (q, p) e R. 
But, recall that h q = v/(c - u), which describes the slope of the vector field in the Eulerian formulation. 
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It follows that, if we restrict C,,i c to be solutions with e > and sufficiently small, then the wave profile, 
indeed all of the streamlines above the bed, are monotonic between crestline and troughline. This implies 
properties (ii)-(iii). □ 

Remark 4.9 Elements of the mixed solution curves C? , andC',, can likewise be characterized 

3,loc 4,loc 

h(q, p) = H(p; A) + e 0i (p) cos q + ^foip) cos (n 2 q)) + o(e), in X, 
where f i, & e R x are found from ( |4.37[ ). □ 



5 Global bifurcation from simple eigenvalues 



The next step in our program is to continue the local curves of Theorem 3.1 and Theorem 1 1 .4| globally. 



In this section we shall address the former (i.e., the simple bifurcation case), while the latter will be the 
subject of section [6] 

As in previous section, we denote 

X := [h e C^(fi) : h = on B), Y = Y x x Y 2 := C£?(R) x C£ T a (T), 

where the subscript "per" indicates 27r-periodicity and evenness in q, R is the rectangle (0, 2n) X (po, 0) and 
T - [0, 2n] X [p — 0}. In this section, we let the nonlinear operator Q — (@\ , Q 2 ) ■ R X X — » Y be given by 

Q x {h) := (1 + hl)h pp + Kohl - 2h g h p h pq - gQi - d(h)) Pp hl + hlB(-p) (5.1) 



2 (Q, h) ■= (l + h 2 q + h 2 p (2o-K[h] + 2gph - 0) 



p=0 



(5.2) 



Recall that k denotes the mean curvature: 

h„ 



K[h] := 



Note that, by construction, the laminar flow solutions //(•; Q) found in Lemma [372] satisfy 

0(H(-; A), Q(A))) = 0, for all A > -2B min + e . 

The main engine for our continuation argument will be Kielhofer degree theory, a generalization of the 
classical theory of Leray-Schauder (cf. 0161 ). The basic outline for this type of argument is by now well 
established (see, e.g., Il6lll4ll27ll36ll ). Fixing 8 > 0, put 

O s := {(g, h) e R x X : h p >6 in R, Q- 2aK[h] - 2gph > 6 

Likewise, denote 

S s := closure in R x X of {(Q,h) e O s : 0(Q,h) = 0, h q $ 

and let C' g be the component of Ss that contains the point (Q*,//„), where 2» := Q(/!„), //, := H(-;A*) 
is the point where Crandall-Rabinowitz bifurcation occurs. We note that this is valid since the laminar 
flows all have zero curvature, hence (Q*,H*) e n^Os- Thus C' s contains the local curve C[ . A typical 
bifurcation diagram for this scenario is shown in Figure 5.1 Note that, the continuum C' s is connected, but 
need not be path-connected in general. We shall discuss this important fact more in the next section. 
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" X 

Os 




Fig. 5.1 

The main result of this section is the following. 

Theorem 5.1 (Global Bifurcation from Simple Eigenvalues) Let 6 > be given. One of the following 
alternatives must hold. 

(i) C' s is unbounded in R X X. 

(ii) C' 6 contains another trivial point (Q, H) e T , with Q + Q*. 

(iii) C' 6 contains a point ( Q, h) e 80s- 

The majority of the work in proving THEOREM |5.1| is showing the admissibility of Q in the sense of Kielhofer 
degree. This will, of course, hinge on being able to derive suitable a priori estimates for the linearized 
operators — a task that is difficult in two respects. First, the addition of the nonlocal operator d to the 
interior equation renders (at least initially) the usual elliptic regularity results inapplicable. The method 
for resolving this problem, first presented H361 . is to "freeze" d to a constant, apply regularity theory to the 
resulting elliptic problem, then argue back to estimates on the unfrozen problem. 

This leads us to the second difficulty: because of the higher order (tangential) derivatives on T, even 
the frozen elliptic problem lies outside the domain of standard Schauder theory. Though slightly more 
obscure, this type of boundary condition is usually referred to as Venttsel-type, and has received some 
attention in the literature (see, e.g. |2T1|221 ). In the interest of self-containedness, we present here one of 
the more basic linear result, as that is all that we shall require in this section. This is done by fleshing out 
a remark in [ 19] on the solvability of such problems in the constant coefficient case. 

Lemma 5.2 Fix n>2 and let Lbe a second-order, linear, uniformly elliptic differential operator on R", 

n n 

L := ^ a u (x)d Xi d Xj + Y bi(x)d Xl + c(x), x e R^; 
i,j=\ i=\ 

M a homogeneous, second-order, linear, uniformly elliptic differential operator on <5R", 

n-l 
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and N a first-order, linear, oblique differential operator on <9R", 

n 

iV := Y Hx')d Xi , x e <5R" , |b„| > j3 > 0. 
i=i 

Suppose that 

lk'./llc 2 +"(K'> c|| c «(Kj), \\o-ij, hWcmi) < K 

and let k^, km denote the ellipticity constants ofL and M respectively. Then there exists C > (depending 
on K,/3, n, Kl and Km ) such that, for all u e C 2+ff (R" ), 

IMIc^dj) ^ C (llwllc^ + \\Lu\\ ca(Wl) + ||(M + N)u\\ , (dK) ) . (5.3) 



Proof First note that, without loss of generality, we may take b-„ c = 0. To see this, first suppose we had 
the lemma for homogeneous L. Then, for general L, we can treat the lower-order terms as a right-hand 
side and apply \53\ to find, 

i 

+ llc|lc<.(Rj)ll"llc«TO + W( M + N)U\\ C ^8RI)\ 

Next, interpolate the C 1+ff -norm on the right-hand side between C a and C 2+a and absorb the latter term 
into the left-hand side to get 

IMIc^CRj) < ^\\u\\ c2+ , r(M) + C (iMI^Rj) + \\Lu\\ cam + ||(M + N)u\\c«(0K>) ■ 
< C(||i/|| c „ (i j ;i + \\Lu\\ or(M) + ||(M + N)u\\ c «(dK)) > 

where, as required, C = C(K,f3, n, Kl, Km). So, indeed, it suffices to assume that L is homogeneous. 

Now let us introduce a few notational conventions. Let the Poisson map, f : C"(R") X C 2+Q, (<9R") — > 
C 2+Q, (R^), be defined by 9(f,g) := m, where m is the unique solution of the Dirichlet problem 

(1 - L)u = /, in R", u(x', 0) = g(x'), Vx' e 3R+. 

By standard Schauder theory this operator is well-defined, bounded and linear in each component. We 
may therefore define a Dirichlet-to-Neumann operator, N, by 

WS):=(3,,WS)) 

Claim. There exists a positive constant C = C(#, n, k l ) such that, for any (/, g) e COR?) x C 1+a (3R" ), 

l|jv(/,s)llc« ( i») < cfii/ii^ + yic.* (flR .,). 
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Proof of Claim. Let g e C 2+a (dR'D = C^R"" 1 ) and / e C a (R" + ) be given. By linearity of L, N(J,g) = 
N(f,0) + N(0,g). The first term, however, can be easily controlled since P(/, 0) satisfies the classical 
Schauder estimate for vanishing boundary data: 

IW.0)H c2 « (i r ) <C||/|| co( ^ )) (5.4) 
where C = C(||ay||c«, «z,,«) > 0. Therefore, 

IIAr(/,g)||c«c»j) < ll//(o,*)||c-(flR-) + \\N(fM\c-mi) 

<\\m,g)\\c° W + \\P(fM\ c ^m) 
<c(||Ar(0,g)l| c . (flR . ) + 11/11^). 

Our main task is thus to estimate N(0,g). Denote by A' the Laplacian operator on R"~'. Suppose at 
first that the coefficients of L are constant. The advantage of this is that it implies 1 —L and 1 —A' commute, 
from which we obtain the following identity 

P(0,(l-/)g) = (l-/)P(0,g). 

Hence, 

||?> (0,(1 -A')g)\\ cUam < ||1 - ^o^o^WM^m 
< C||g||c3+«(«y 

<C\\{\-A')g\\ 0+a(mi) , (5.5) 

where the constant C - CQ\aij\\c,KL,n) > 0. 

Standard elliptic theory tells us that 1 - A' is a bounded linear operator mapping the space C 2+a (R" _1 ) 
surjectively to C Qr (R"~ 1 ). Then, replacing g by (1 - A')~ l g in ( |5.5| l, we get 

IP(0,g)ll c , t . S) < Cyici^flR.). (5.6) 

The claim is an immediate consequence of this inequality, since 

\\N(o,g)\\ c « mi) < \\no,g)\\ c ^ (W) < c\\g\y +a(aK) . 

Combing this with ( |5.4[ l we have the desired result in the constant coefficient case. 

In order to complete the proof of the claim, therefore, we need only to reprove ( |5.6[ ) for variable 
coefficients. Naturally, the idea will be to estimate the commutator [1 — L, 1 - A'] effectively. 

Put u := f{0, (1 - A')g), v := ^(0,^). That is, unraveling notation, take u and v to satisfy: 

(1 - L)u = 0, in R", w = (1 — A')g, on <9R" , (5.7) 

(1 - L)v = 0, in R" , v = g, on 31^. (5.8) 
Put w :- (1 - zl')v. Then, from ( |5.8| l, we have that w satisfies 

(1 - L)w = [1 - L, 1 - A']v, in R" , w = (1 -A')g, on 3R^. (5.9) 

Here we have let [■, •] denote the commutator. Subtracting \5.1\ from \5.9\ leads to 

(1 - L)(w - u) = [1 - L, 1 - zl']v, in R^, w - u = 0, on 3R^. 
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Applying the Schauder estimates to this last equation we find 

\\w - u\\ c2+am < Cl|[l - L, 1 -^'Ivllc.^. (5.10) 

Here C = C(||ay||c*,<L»«) > 0. But note that [1 — L, 1 - zf] is a third-order operator whose coefficients 
involves up to second-order derivatives of a,j. We therefore obtain from estimate ( |5.10| l that 

||ii - (1 - 4')v\\ c2+llWl) < C\\v\y + „ m) < C\\g\\ C 3» idK) < C\\(\ - A')g\\ c i» mi) , (5.11) 

where the constant C = C(K, k l , n) > 0. 

Now, rephrasing this in the earlier notation, we find 



||no, (i-^)|| cl ^< lla-Wo, 



+ ||y(0, (1 -A')g) - (1 - A')P(0,g)\\ cl+am . 



The first term on the right-hand side can be estimated as in ( |5.5| ), while the second term is bounded from 
above by ( |5.1 l| l. Equivalently, if we replace g by (1 - A')~ l g, it follows that 

\W(0,g)\\ o , {aK) < \\P(.0,g)\\ c u« m) < C\\g\\ c n a(dK , 

where C = C(K, k l , n) > 0. This proves the claim. □ 
Let u e C 2+a (R") be given. For notational convenience, we let h denotes the trace of u on 3R". The 
standard Schauder estimates for the Dirichlet problem are 

IMIc^R?) ^ c (iMIcacj + PIIc*~wri> + ll L "llc««)) ' (5 - 12) 

where the constant C - C(\\aij\\ c «,KL,n). 

Let (/, g) e C 2+a (Wl) X C a (dM" + ) be given, and suppose 

(1 -L)u =f, inR" + , (M + N)u = g, ondR n + . 

We can express the second equation above entirely as a problem on 3R" by means of the operator N( f, ■): 



' n-\ n-1 ^ 

\i,j=l i=l 



h = -b„N(f, h) + g, on 3R" . (5.13) 



Suppose that n > 3. Applying the Schauder estimates on R" 1 to this equation yields 

Plb«(*o < C (\\h\\ c « {dK) + \\N(f, h)\\ C r (aK) + Motin®) , (5-14) 

where C - C(K,K M ,j3,n) > 0. Applying the claim to the TV-term above and interpolating the resulting 
C 1+ff -norm between C a and C 2+a we obtain 

Pllc 2 +"(«K[) < C (\\h\\ <(dKl) + ll/ll C «(R») + Plb+°(«;) + llgllc«C3RJ)) 

< ^\\h\\ c ^(3w + ) + c (WMctdEZ) + H/llc»«) + HsIlccflRJ)) 

< C (\\h\\ c „ (dK) + \\f\\ c „ (Wi) + \\g\\c« m v) ■ (5.15) 

Here C = C(K,K L ,K M ,/3,n). 
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Inserting ( |5.15[ ) into ( |5.12| l — keeping in mind that u = h on <9R" and ll/l^mr-, < \\u \\ c „^ + \\Lu\\ ca ^ 
— we arrive finally at the following: 

IMb««) < c(||u|| ca(S 5- ) + \\Lu\\ cam + \\(M + N)u\\c-(dK)) • ( 5 - 16 ) 

where C = C(A", k m ,/3, n). This completes the lemma for n > 2. 

For the case n — 2 we introduce an artificial variable, ?, to move up a dimension. Notice that the 
only place where the previous argument fails is in deriving ( |5.14[ ), as Schauder theory would require 
the boundary to be at least two-dimensional. Let h, f, g and N( f, li) be as above and let their extensions 
(h,g,N(f,h)) € C 2+£ *(R 2 ) x C a (R 2 ) x C 1+a (R 2 ) be defined as follows 



h(t, x) := h(x), g(t, x) := g(x), N(f, h){t, x) := N(f, h)(x), V(t, x) e R z . 

Of course, for any k > 0, the C* +ff (R 2 )-norm of the extensions (h, g, N) coincides with the C Ar+a (R)-norm 
of the original original functions (h, g, N). 
Define operators M, by 

M:=M + d z , N = N + d t . 
Then (M + N)g - h iff (M + N)g = h, so from ( |5.13| l we obtain 

{a u d 2 . + d 2 + b t d x + d t ) h = -b 2 N(f,h) + g, on R 2 . 

As we are now in the proper dimension, the usual elliptic estimates yield ( |5.14| l with h, g and N in place 
of h, g and N, respectively. Then, using the fact that the extensions have the same relevant norms as the 
original functions, we may argue as before to get ( |5.15[ ) and thereby reestablish ( |5.16| l for n — 2. □ 



Remark 5.3 A standard cut-off function argument generalizes Lemma [5T2| to the case at hand, i.e. where 
we have Q — R, is bounded with periodic boundary conditions on the lateral sides and Dirichlet data on 
the top. □ 

Since we now have adequate a priori estimates for the frozen problem, we fix any t e R and define 
@ {T) : R x X -> Y by 

g[ T \h) := (1 + h 2 q )h pp + h qq h\ - 2h q h p h pq - g(h - T)h 3 pPp + hl/3(-p) (5.17) 



6%\Q,h) := & 2 {Q,h) = (l +h 2 q + h 2 p (2<r«[h] + 2gph - g)) 



(5.18) 

p=0 



We will show that, for h £ Og, is a uniformly elliptic, quasilinear differential operator, while {p? ls 
a non-degenerate, uniformly oblique Venttsel-type boundary operator. For later reference we now compute 
the Frechet derivatives of Q and Q ij) : 

lh (h) = 2h q h pp d q + (1 + h 2 q )d 2 p + 2h p h qq d„ + h 2 d 2 q 

- 2(h pq h q d p + h pq h p d q + h p h q d p d q ) (5.19) 
- 3g Pp (h - d{h))h 2 p d p - g Pp h 3 p (\ -d) + 3h 2 p(-p)d p 



02h(Q,h) = (2V, + 2h p (2a«[h] + 2gph - Q)d p + 2gph\ + 2o-h 2 /[h]j 



(5.20) 
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&f h (h) = g lh (h) + 3 g (d(h) - oytjPpdp + gp P K d 



= 2h q h pp d q + (1 + h L q )d p + 2h p h qq d p + h 2 p d 2 q 
- 2{h pq h q d p + h pq h p d q + h p h q d p d q ) 

- 3gP P (h - T)h 2 p d p - gp p h 3 p + 3h 2 pP(-p)d p 

6 { :l(Q,h)=&2h(Q,h), 



(5.21) 



(5.22) 



where 



k'[K\ := 3h qq (l + h 2 q \ 512 d q - (l + tf q y £ d 2 . 
, we now establish the key admissibility lemmas. 



-3/2 



Following 



Lemma 5.4 (Proper Map) Suppose K is a compact subset ofY and D is a closed, bounded set in Og, then 
Q~ l {K) n D is compact in RxX. 

Proof Let {(fj,gj)} be a convergent sequence in Y = Y\ X Y2 and assume that for each j > 1, (fj,gj) = 

@{Qj,hj) for some (Qj,hj) e O s with [hj] bounded in C^(R), Q, bounded in R. We wish to show that 
there exists a subsequence of {(Q j, hj)} convergent in R X X. 

Denote 8j := d q hj, for j > 1. Then, differentiating the relation between (Qj, hj) and (fj, gj), we find by 
( |5TT| > that, for all j>\, 

d q fj = d q g x {hj) 

= (1 + (d q hj) 2 )d 2 p 6j + (d q (l + (d q hj) 2 ))d 2 p hj + (d p hj) 2 d 2 0j + (d 2 q hj)d q (d p hj) 2 

- 2(d q hj)(d p hj)d p d q ej - 2(d 2 q hj)(d p hj)(d p d q hj) - 2(d q hj)(d p d q hj) 2 

- g(d p hj) 3 Pp 8j - g(hj - d(hj))d q (d p hj) 3 + d q (d p hjf/3(-p), 
Grouping terms above we may rewrite this in the form 

(l + {d 2 q hj) 2 )d 2 p ej + (d p hj) 2 d 2 q 0j - 2(d q hj)(d p hj)d q d p 0j = 

d q fj + F(d q hj, d p d q hj, d 2 hj, d p hj, d 2 p hj, d(h)), in R, 



in R. 



(5.23) 



where F is the cubic polynomial dictated by the previous expression. Notice that the coefficients of the 
higher order terms above satisfy: 



4(1+^2-4^ 



4h 2 p > S 2 . 



Hence, the right-hand side can be viewed as a uniformly elliptic operator acting on 8. 
We may likewise differentiate the equation on T to discover that, for each j > 1, 

dqgj = d q @ 2 {Qj,hj) 

= 2(d q hj)d q 8j + 2(d p hj)(2o-K[hj] + 2gphj - Qj)d p 0j 
+ 2gp(d p hj) 2 h q + 2o-(d p hj) 2 (d q K[hj]), on T. 

+ 2(d q hj)d q 0j + 2(d p hj)(2o-K[hj] + 2gphj - Qj)d p 0j = 

d q gj +G(d q hj,d p hj), 



Equivalently, 



2(T(d p hj) 2 d q 



d q 6j 



(5.24) 



on T, 
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where G is the quadratic polynomial determined by the previous equation. Observe that the boundary 
operator present on the left-hand side above consists of a one-dimensional elliptic operator added to a 
uniformly oblique operator, since 

\(2o-k[H] + 2gph - Q)h p \ >S 2 on T 
for (Q, h) e O5. Finally, we note that 

9 j = on B. (5.25) 

Since {(/}, gj)} is convergent in C^"(R) X C^"(T) we have that the sequences {d q fj}, {d q g j] are Cauchy 
in Cp er (R) and Cp er (r) respectively. Now, by assumption, {hj} is uniformly bounded in C^"(R). From this 
it follows that F(d q hj,d p d q hj,d 2 hj,d p hj,d 2 hj,d(h)) is uniformly bounded in f(R), and G(d q hj,d p hj) 
is uniformly bounded in C 2 ^"{T). Each of these, viewed as a sequence in j, is therefore Cauchy. Then, 
by the compactness of the embeddings, the right-hand sides of equations ( |5.23[ > and ( |5 .24[ > are strongly 
pre-compact in Cp ei .(R) and Cp er (r), respectively. Possibly passing to a subsequence, we may take both to 
be convergent in these spaces. 

We now consider differences 6j — Ok, for j, k > 1. By ( |5.23[ ) we have 

Fjk = (1 + (d 2 q hjf)d 2 p (0j - 6 k ) + (d p hjfd 2 q (0j - k ) 
- 2(d q hj)(d p hj)d q d p (0j - 6 k ), in R, 

where by our arguments in the previous paragraph we know F jk -> in C a (R). Similarly, from ( |5T25] > have 
that 6j - 6k vanishes on the bottom and on the top ( |5.24[ ) tells us 



G jk = 2cr(d p hj) 2 d q 



MOj-Ok)) 



+ 2(d q hj)d q (0j - 6 k ) 



(1 +/j2)3/2 

+ 2(d p hj)(2cTK[hj] - Qj)d p (9j - 9 k ), on T. 

Here Gjk — » in C a (T), again by the considerations of the preceding paragraph. We now apply the mixed- 
boundary condition Schauder estimates of Lemma 5.2 to the differences 9j - 9 k to deduce that there exists 
some constant C > with 

c\\6j - Okiy^ < \\F jk \\ ca(1) + WGflMcm + llfj - Wcm- 

Thus, Qj - 9 h — > in C 2+a as j, k — » 00 . We have shown, therefore, that all third derivatives of {hj} are 
Cauchy in Cp er (R), except possibly for [dphj}. To demonstrate the same holds for {d p hj}, we use the PDE 
to express d 2 hj in terms of the derivatives of order less than or equal to two: 

d 2 p hj = (1 + (d q hjfr l (fj - d 2 q hj{d p h) 2 + 2{d q h){d p hj){d p d q hj) 

- g(hj - d{hj))p p {d p hjf + (d p hj) 3 B(-p)), in R. 

But we have seen that the right-hand side is Cauchy in C 1+ °(R), hence [dphj] is also Cauchy in Cp er (R). 
We conclude that the original sequence, {hj}, is Cauchy in C 3+ff . □ 

Lemma 5.5 (Fredholm Map) For each (Q, h) e Og, the linearized operator @h(Q, h) is a Fredholm map 
of index Ofrom X to Y. 



48 



Samuel Walsh 



Proof Let if/ e C^ a (R) be given and fix t, QeR,heX. Put 

(O := d q (gff(Q, h)ffl) - (d q &%(Q, h)) [<fr], for i = 1,2. 

Here, by (d q Q ( ^(Q, h))[if/] we mean differentiating the coefficients of Q^) in q, then applying the resulting 
operator to if/. Then d q if/ satisfies: 

Q^(h)d q if/ = 0« on/?, 
g%(jQ,h)d q ilr = 0® onT, 

c^iA = on 6, 

which is a uniformly elliptic PDE with a Venttsel boundary condition. Thus there exists a constant C > 0, 
independent of t/f, such that 

<WHc*«<K) < lia^llc.^ + II0 (1) |I C "(«) + U (2) \\0'(T) 



< 



ll C 2 + «(«) + lia ? (^u(A¥)llc»® + K(^2(e,/#)llc«cr) 



On the other hand, we may express c^i/f via the partial differential equation to arrive at an estimate for dpif/ 
of the same type. Combining these estimates we find that, for some C > and all if/ e X, 



C\m\c^m * \Wc^ ( R) + W^Wcw + W (2) Wc<(T)- (5.26) 



If we now apply Lemma 5.2 to if/ directly, we can estimate the C 2+ff -norm of if/. Then ( |5.26| l becomes 

CUWx < UWy, + llMfcWlk + II0»(G. Willi- (5-27) 

This is the key estimate, from which we can conclude the range of @r?? is closed and the null space 
finite dimensional. The task at hand, therefore, is to prove it with Q in place of Q (T) . As in ll36ll . this is 
made possible by the following trivial but convenient fact: for any if/ as above we have 



\dffl\ < £ \if/\dq < UWc^R)- 



We are therefore able to estimate terms involving d easily in spaces of Holder continuous functions. In 
particular, 



\\g lh (h)if/ - ffS^^cm = W^i.d(h) - rWpPpif/p + gp P h p dffl\\ ca(1) 

< Q\if/\\ cl+am , (5.28) 

where the constant C above depends only on p p , \\h\\x and our choice of cr. Likewise, an identical argument 
gives the estimates 

\m u (h)if/ - di@VQi)if/\\ cail) < C\\if/\\ c2+ „ (1) , i = p,q 
where C depends again on p p , \\h\\x and the choice of t. Of course we do not need to make such arguments 
to estimate the (f^Qi) term, as it identical to @2h(h). 

Combining these observations with ( |5.27[ l, we find that, for some C > and all if/ e X, 

CU\\x < M\ci~m + WSnOm* + \W2h(Q,h)if/\\r 2 - (5.29) 
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Applying the arguments of 1361 . we are able to conclude that @i,(Q, h) is Fredholm. The Fredholm index 
has a discrete range, hence by the connectedness of 0$ it must be constant on this set. But @h(Q*,H„) 
was shown to have a one-dimensional null space and range of codimension one. Since {Q,, H t ) e Og by 
construction, the Fredholm index is uniformly along the continuum C' s . □ 

Lemma 5.6 (Spectral Properties) (i) V6, e > 0, 3c\ , c 2 > such that for all (Q, h) e Os with \Q\ + \\h\\x < 
M, for allifj € X and for all p e C with \p\ > c 2 and |arg /x\ < n/2 - e, we have 

CiUWx < M o/2 1ICA-Ai)0l7, +l//r /2 ||(B-^lly 2 
where A = A(Q, h) = Q xh (h) and B = B(Q, h) = 2h {Q, h). 

(ii) Let E - E(Q, h) denote the spectrum of (A, B). Then X consists entirely of the eigenvalues of finite 
multiplicity with no finite accumulation points. Furthermore, there is a neighborhood N of[0, +00) in the 
complex plane such that Z (A, w) n N is a finite set. 

Proof To prove (i) we follow the argument of Agmon (cf. HI). The Venttsel boundary conditions, however, 
necessitate some modification of the standard argument. All that we must do here is trivially modify the 
technical tricks of that argument. Part (ii) follows immediately from part (i) in light of the previous lemmas, 
so we omit the details. 

Fix r e R with |r| < M, and let tfr e X, p € C be given. Put 8 := argp, and suppose for some e > 0, 
\8\ < n/2 + e. Consider the operators 

A (T) := £ ( ^(/z), Df := A {T) + e w cr}, on R x R, 

and 

Df = D 2 := B + e ie d], on R X T. 

Then, for any t, is elliptic with constant of ellipticity independent of t, while is of the form of 
Lemma 5.2 Let ( : R — » R be a cutoff function supported compactly in the interval / := (-1, 1). Put 

e{t) := e^' Zt m <f>(t,q,p) ■= e(t)i/,(q, p). 



We apply the Schauder estimates of Lemma 5.2 in R 3 with the boundary operator Dj to <p{t, q, p) to deduce 



2 



C||0lb +W < \W\o<(ixR) + Wi'HcexR) + ll£>2 >Hc«(/xr) (5.30) 
for some constant C > inde pend ent of ifj. In fact, since t occurs as the coefficient of a first-order term in 



&il(Q> h), we have by Lemma 5.2 that C can be chosen independently of r (recall, |t| < M.) 



A quick calculation readily confirms the existence of C, C > 0, depending only on with 

Cp a/2 < ||e||c«(R) < C'p al \ Cp {X+a)l2 < \\e\\ c u„ m < C'p (l+a)/2 . (5.31) 
Using this, we can unpack ( |5.30| l to derive the following set of estimates: 

< c\p\ al2 - m\c(fD + mihW - eftmncw) 

+ ii(/N 1/ V+r)A Ml/2 vii CTO 

< C'W 12 (\\(0 lh (h) - p) n c <( R ) + IMIc»cr>) (5-32) 

+ C'(^| (1+a)/2 + |//r /2 )|| ( Allc«w (5.33) 
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for \fi\ sufficiently large. In transitioning from ( 5.32[ ) to ( |5.33[ l, we have made use of ( |5.28| l taking t = d(h) 
(which is valid since d(h) < \\h\\x < M and ( 5.30| > holds for any t.) A similar analysis of the boundary 
terms yields, 



C<*(/xT) 



< C'M a/2 \\ (g 2 h(Q, h) - fi) M\c<n + C (l^l (1+£?)/2 + W 2 ) UWoht), (5.34) 

for is sufficiently large. 
On the other hand, 

2,a;RxI — ,a\Rxl + [d t e(t)d q l//]o :a -Rxl + [d t e(t)d p l//]o,a;RxI 

,a;Rxl + 

[e(t)d 2 p ifr] Q 

,a;RxI 

> C"\ia\[l//]o, a ; R + [fih,a-Jt, 

UWCHIXR) > C"\fl\\\n C 0(R) + \M\c\R)> ( 5 - 35 ) 

for \fi\ sufficiently large, and for some C" > independent of /i and tf/. Here, [-]k, a jxR denotes the usual 
Holder semi-norm. Combining ( |5.30| ), ( |5.33| ), ( |5.34| l, and ( |5.35| l, we find derive the following: there exists 
c\,C2 > such that, 

ciU\\^(R) * \A al2 W (A - M\c® + l-"r /2 H ( B - l*M\cm ' 

for all /i e C with > c 2 . If we now repeat the same argument with iff q and if/ p in place of iff, we arrive at 
the inequality in (i). □ 



In order to complete the proof of Theorem 5.1 we will make use of a generalization of Leray-Schauder 
degree due to Kielofer (cf. 01611171 ). With that in mind, we recall the following definitions. 

Definition 5.7 (See Definition II.5.2 of ifTTIl .) Let W, Z be real Banach spaces with W continuously em- 
bedded in Z and let 'W be an open, bounded subset of W. We say that a map G e C 2 (*W; Y) is admissible 
provided that it satisfies the following: 

(i) G is a proper map. 

(ii) For all w e r W, the linearized operator G w (w) is Fredholm of index 0. 

(iii) There exists e e (0,tt/2), a e (0, 1), and constants c\,c 2 > such that, for all ift e X, w e f W, and 
fi e C with |ju] > C2, |arg ii\ < n/2 - e, we have 

cill<Alk<N ff/2 ||(G lv (w)-^|| z , 

where we have complexified W and Z. 

(iv) For each w e f W, the spectrum E of G w (yv) consists entirely of eigenvalues of finite multiplicity with 
no finite accumulation points. Furthermore, there is a neighborhood N of [0, +oo) in the complex plane 
such that Z(A, w) n N is a finite set. 

The significance of admissible maps is that they permit us to define a degree. Suppose that G is such a 
map and let z e Z \ G{&W) be a regular value of G. We define the (Kielofer) degree of G at z with respect 
to < W to be the quantity 

deg(G,1V,z):= £ (-l)" w , 



Water waves with surface tension 



51 



where y(w ) is the number of positive real eigenvalues counted by multiplicity of G w (w). By part (iii) of 
Definition 5.7 we have that v(w) is finite. Likewise, part (i) ensures G~ 1 ({z})Ci'W is a finite set. The degree 
is therefore well-defined at proper values. In the usual way, this definition is extended to critical values via 
the Sard-Smale theorem (this is valid by since G is Fredholm.) 

Our interest in degree theory stems from the fact that the degree is invariant under any homotopy that 
respects the boundaries of IV in the following sense. Let It c [0, 1] X < W be open. Define U t :— {w e 
W : (/, w) e 1/} and likewise dll, := {w e <W : (t, w) e dU\. 

Lemma 5.8 (Homotopy Invariance) The degree is invariant under admissible homotopies. That is, sup- 
pose H c [0, 1] x *W is open. IfH e C (14; Y) is proper and, for each t e [0, 1], fiit, •) is admissible, 
then we say H is an admissible homotopy and have 

deg (7/(0, -),U ,y) = deg (7/(1 , •), *Uu y) 
provided y £ Kit, d!4 t )for all t e [0, 1]. 

The reader is directed to 1171 for a thorough treatment of degree theory and a proof of this result. 
We are now prepared to prove the main theorem of this section. 



Proof of Theorem 5.1 From part (iii) of Definition 5.7 we see that in order to apply the Kielhofer degree 



theory, we must first make a technical change. For h e C^"(R), put 



h:=(h,h\ T )eC£?(K)xC£?(T), 



and 



X := ((Mir) : h e C^(R)} c c£?(*) x C^{T). 
If we let n be the bijective projection n : X — » X, h ^„ (h, h\j), then we can likewise define 



O s := MO s ), 



C' s = n(C' s ). 



Let TV be an open bounded subset of X and fix Q e R. By virtue of Le mma |5 ,4| -Lemma |5 ,6| we conclude 

Y is admissible in the sense of Definition : 



5.7 



This is a trivial consequence 



that the map @(Q, •) : n(-W) 
of the fact that n is a smooth bijection from X onto X. 

From our p revio us analysis it follows that Q(Q\, •) t-> QiQ.2, ) is an admissible homotopy. Thus, in 
light of Lemma 5.8 the degree will remain constant as we move along the continuum C' s . Assuming that 



all three alternatives of Theorem 5.1 fail, we use this feature to generate a contradiction. At this stage, we 
have reestablished all the relevant properties of Q that were true in the constant density case. Repeating 
(verbatim) the proof in J51 — with Q in place of Q — we obtain the theorem. □ 

Remark 5.9 The motivation for choosing the Kielhofer degree (rather than that of Healey and Simpson 



used, for example, in |6 |) lies in Lemma 5.6 due to the presence of higher-order derivatives, the boundary 
operator and interior operator are equivalent, in the sense of the Agmon-type spectral estimates of Lemma 
5.6 (i). This observation was first made by Wahlen in his analysis of the pure capillary wave case, where 
he also introduced the clever trick of considering h, X, etc. (cf. II35I ). □ 
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6 Analytic global bifurcation 

In the previous section we obtained global bifurcation results for the case when the eigenvalue at (//„, Q,) 
is simple. The basis for that method of argument, which harkens back to Rabinowitz 0271 . is the topological 
theory of degree. In the present section, we present an alternative theory of bifurcation that instead relies 
on the (real) analytic properties of the nonlinear operator Q. 

This approach is originally due to Dancer (cf. SHU, or Q for a survey) and has recently been applied 
with great success to the study of Stokes waves (cf. [2 3 4 1). The major insight of Dancer's theory is that 
the germ of the zero-set of an analytic operator has an extremely special structure. Within the context of 
global continuation this means, for example, that if C is an analytic curve that corresponds to the zero-set 
of an analytic function, then there is a unique, maximal, analytic extension of C. Lifting this fact to infinite 
dimensions (and assuming some additional compactness properties) provides a very powerful set of tools 
for global bifurcation theory. 

With respect to the degree theoretic methods of the previous section, Dancer's analytic continuation 
method offers several advantages. The first of these is that it provides a more natural way to treat the 
double bifurcation case. Note that, as the multiplicity of the eigenvalue will be even, a Rabinowitz-type 
argument requires some additional work; in particular, special attention must be paid to the scenario where 
the continuations of the local branches intersect. On the other hand, since maximal analytic extensions are 
unique, when working within the analytic framework, we are untroubled by such intersections. Secondly, 
since we are using more than simply topological information, we obtain considerably better regularity for 
the global bifurcation branch. Indeed, we will prove that the solution continuum is, in fact, path-connected 
and admits a locally injective, analytic reparameterization (in the sense of Definition [TT|) 

Let us begin by presenting the central result of the analytic theory of global bifurcation. This was 
originally proved by Dancer (cf. fi WTTIO . but we shall use a slight rephrasing of that result due to Buffoni 
and Toland (cf. 0). For the purposes of this exposition, we shall temporarily set aside our previous 
notation. Let real Banach spaces W, Z, open set U C R X W and a real analytic T : U — > Z be given. 
(Looking forward, we should think of U as taking the place of 0$, W that of X, Z that of Y and T that of 
Q.) If the following structural properties hold 

(Gl) T := {(Q, 0) : Q e R} c U and T(Q, 0) = 0, Vg e R, 

(G2) f w (Q, w) is a Fredholm operator of index zero for each (Q, w) e U such that T{Q, w) = 0, 
(G3) 3g, e R such that (i) N(T w (Q*,0)) = [sfa : s e R), and (ii) T^Q,, 0)0* t K{r w {Q„ 0)), 



then, by the analytic equivalent of Theorem |3. 1 1 there exists an analytic local bifurcation curve at (<2„, 0). 



That is, there exists an e > and an analytic (Q, w) : (-e, e) — > R X W such that !F(Q(s), xd(s)) = 0, for all 
5 e (-e, e), Q(0) = Q* and vo'(s) = (f>*. Let us denote 

C; 0C :={(Zi(s),w(s)): se(0,e)}, S := {(Q,w) e U : T(Q,w) = 0). (6.1) 

By taking e sufficiently small, we may assume that m' + and C* c S\7~. 

Theorem 6.1 (Analytic Global Simple Bifurcation) Let W, Z, U and T be given as above and assume 
that (Gl), (G2) and (G3) hold. We may then take (Q, m) as before and let C^ oc and S be defined by \6.\\ . 
Suppose further that 

Q' i Oon (-£,£), (G4) 
all bounded closed subsets of S are compact. (G5) 
Then there exists a unique, maximal, continuous curve "K that extends Cj 1 " in the following sense. 
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(a) 7C = {(CO), to(s)) : 5 e [0, oo)) c I/, w/zere (G, ro) : [0,m)^1xW is continuous. 

(b) C+ c c 7C c S. 

(c) 77?e wf {s > : A^(^,(JQ(s), to(s))) # {0}} /icw no accumulation points. 

(d) A? eac/z point, 'K has a local analytic reparameterization. 

(e) One of the following occurs. 

(i) ||(jQ(j), V0(s)\\ ->oofljS-> oo; or 

(ii) (fi(.s), tD(i)) approaches the boundary ofUass^> oo; or 

(iii) 7C « a closed loop, i.e. 3T > .smc/j f/zaf + T), trj(,s + T)) = (Q(s), w(s)), for all s e [0, oo). 
This implies that 7C = {(£>(*), to(») : 5 e [0, T]}. 

(f) If, for some < *i < S2> wc Ziave 

(JQCti),»Cti)) = (G(52),tD(j 2 )), anrfAr(Tw(JD(Ji), »(«))) = (0), 
then (e)(iii) occurs and S2 - s\ is an integer multiple ofT. 

In order to apply this theory to our problem we need to make several benign alterations to hypotheses 
(G1)-(G5). First, we wish to replace T occurring in (Gl) with that from the previous sections, namely the 
1 -parameter family of laminar flows. This is justified since the particular form of the trivial solutions plays 



no role in the proof of Theorem 6. 1 Indeed, once the local curve Cj 1 " is found, no more mention need be 
made of T. 

Secondly, we wish to consider the case where double bifurcation occurs at the local level. Again, this 
is possible because the purpose of (G3) is merely to guarantee the existence of C^ c and has no bearing on 
the global analysis. If, as in section|4j there is a set of local bifurcation curves {C\ocj}^ =v then there exists 
a corresponding set of local parametrizations {(&;, H);)K=i anc ^ s0 we ma y define {CjL ■}■_ l in the obvious 



way. Applying the argument of Theorem 6.1 to each of these curves yields the following corollary. 



Corollary 6.2 (Analytic Global Double Bifurcation) Let W,Z, U and T be given as before. Suppose that 
the following structural properties are satisfied: 

there exists a family of trivial solutions T c U n S, (GT) 

T W (Q, w) is a Fredholm operator of index zero for each (Q, w) e U n S, (G2') 

3(Q t ,w f ) e T, e > 0, and analytic functions {(Q,-, tx>j) : (-e, e)^Rx W}f =1 

such that, for all s e (— e, e), i = 1,. . . ,N, (G3') 

(Q,(0), ro,(0)) = (jg.,w.) andT(£n(s), Wi(s)) = 0, 

Q; £ on (-e, e), i=l,...,N, (G4') 

all bounded closed subset of S are compact. (G5') 

For each i — l,...,N, define Cjt , : = {(£tj(s),Wj(s)): s e (0, e)}. Then there exists a set of continuous 
curves CK,}^ with "7C the unique maximal analytic extension ofC^ ocj in the sense of Theorem 6.1 for 
i'=l,..., a/. — 
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This result follows directly from the proof of Theorem 6. 1 keeping in mind our observation of the 
last several paragraphs. For brevity, we therefore do not provide a proof. Note, however, that this is not 
a new result: the multiple bifurcation case was considered by Dancer in (5J. Corollary 6.2 is merely a 
translation of Dancer's result to the language of Buffoni and Toland. 

With the background results in place, we now return to the notational conventions of section [5] and 
prove our main result. Fix 8 > and define 

O s := {(G, h) e R x X : h p >6mR, Q- cr K [h] - 2gph > s] , 

S := {(g, h) e R x X : Q(Q, h) = 0}, S 6 := closure in R x X of (S n 0, 5 \ T). 

By the local analysis of the previous sections, we may take the non-laminar solution set near (H*, Q*) to 
be comprised of a set of C 1 -curves, [C[ )._ v for some N > 2. We now prove that these curves are, in fact, 
analytic. 

Lemma 6.3 (Analyticity) (a) The maps A t-> H(-;A) e C 3+a ([po,0]) and A t-> Q(A) e R are analytic on 
the interval (-2B m ;„ + ej), °°)- Furthermore, if I is any subinterval of(-2B m i n + €q, oo) that contains A t , then 
the operators T : / XX — > Y and g:RxX^f (defined by p.20[ )-( [X2T| anrf ( |5 l| l-( |5.2| l, respectively) 
are analytic. 

(b) TTzere exwta e > a«d a sef of analytic functions [(J&i(s), t),-(^)) : (— e, e) h* R x X}^ f/iaf 
q 0Cji = {(Q,(i), Wj)) : s e (-6, e)},/or i = 1, . . . ,AT. 

Proof The analyticity of £r is clear, since 

d j Q Q, die = 0, V; > 2, k > 5. 

Notice also that, since 

0(A) = ^ + 2gp(0)H(0; A), T(A, w) = §{Q(A), w + H(; A)), 

it suffices to prove that A t-> //(•; /l) is analytic in order to obtain part (a). 

We shall do this by extending A i-> //(•; A) to a function defined on the complex half -plane H := \z e 
C : Rez > -2S m i n + eo)> an d then showing that the extended function is complex analytic. 

Mimicking the proof of Lemma 3.1 we can prove that, for each A e H and p e [po, 0], there exists a 
unique H = H(p; A) e C satisfying the following implicit relationship 



H(p;A)= f - — , pe[p .0], AgH, (6.2) 



, P0 + G(s;/l)' 
where we are taking the positive branch of the square root and 



G(p; A) = 2B(p) + 2g f Pp {s) (H(s; A) - H(0; A)) ds. 

J Pa 



If we denote A - A\ + 1A2 and fix p e [po, 0], then 

d lid d \ 

H(p;A)= -\—+i—)H(p;A) 
dA 2\dAi dA 2 J 



= - fig [ p p (j)4: i) - ff(0; i)) rfj) (A + G(r; i))- 3/2 dr. 
"Jpo \ Jpo dA } 
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The last line above was found by differentiating ( |6.2| l (with respect to A\ and A2), and making use of the 
definition of G. 

It follows that, if we denote Y(p; A) :- H(p; A) - H(0; A), then 



d_ 
dA 



Y(p;A)= [ Ig [ p p (s)-^Y(s)ds)(A + G(r,A)r 3/2 dr. 

■Jp \ -J pa dA I 



(6.3) 



Recalling the definition of eo in ( |1.16| l, we see that, since A e '. 



\A + G(p;AT V2 < |Ai + 2S min r 3/2 < (2g\p Q \ 2 \\pp 



,-1 



Inserting this last expression into (|6.3|l, we find 





1 


d 


tYi,A) 


< - 


-=Y(-;A) 


dA 


co " 2 


dA 



and thus dY/dA = 0. As H(p; A) = Y(p; A) - Y(p ; A), this implies further that dH/dA = 0, hence A t-> 
H(-; A) is complex analytic. 

Now, (b) follows directly from the fact that T is analytic. Indeed, simply substituting the Crandall- 
Rabinowitz and Lyapunov-Schmidt theory of the previous sections with their analytic counterparts, we 
obtain the analytic parameterizations {(Q,, l);)}^ as claimed. □ 



Following the notation of Theorem 6.1 let us denote 



C' :=m(s)Ms)):s€(0,e)}, i = 1, 



,N. 



Note that, to make our notation more concise, we have dropped the + superscript. Our main result for this 
section is the following. 

Theorem 6.4 (Analytic Global Bifurcation) Fix 6 > 0. There exists a maximal, global, continuous exten- 
sion 'Kg. c S s ofC' Si (in the sense that Theorem 6.1 a)-(f) holds for < K' Si ), for i — 1, . . . , N. Furthermore, 
for each i, one of the following alternatives must hold. 

(i) ||(AO), I);(s))||rxx -* °° as s -» 00; or 

(ii) the closure ofK' gi contains a point (Q, h) e 80s; or 

(iii) 7C ; is a closed loop. 



Remark 6.5 In several respects, the alternatives above are considerably stronger than their analogues in 
Theorem |5.1| If, recalling the notation of section |5j we define C' g to be the connected component of S s 
containing 0;C[ ., then clearly {Jj'K.g c Cs- In other words, the conclusions of Theorem 6.4 applies 



only to a (very special) subset of C' g , whereas the alternatives of Theorem 5.1 apply to the entirety of C' s . 
In that sense, the theorem above is stronger, since it identifies — in fact, it constructs — the precise path 
along which the alternative occurs. Moreover, the fact that [jj'KL is path-connected and locally analytic 
is noteworthy in and of itself. This cannot be obtained through purely degree theoretic methods, since it 
requires more than topological information about the solution set. 



Still, it is not possible to completely reconcile the alternatives of Theorem 5.1 and Theorem 6.4 For 
instance, it is entirely permissible for ||(Q,(i), h,(i))|| —* 00 as J —* 00, but, for some s > 0, we have 
(£ii(s),\)i(s)) e T. Moreover, it may be that C' Si is a closed loop for all i - l,...,N, yet U;^, n 7" = 
(Q t .,H t ). This is possible, again, because |J, 'KU is only a subset of C' g . 

Finally, let us note that a similar result can be proved for the simple bifurcation case — we need only 
use Theorem |6 . 1 1 directly, rather than Corollary |6.2| □ 
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Proof of Theorem 6.4 In light of Corollary 6.2 and Lemma 6.3 the structure of our argument is fairly 



straightforward: we need only confirm that Q satisfies structural properties (G1')-(G5')- 

With that in mind, let X, Y, 0$ and Q stand in for W, Z, U and T of Corollary 6.2 respectively. Then 



the family of laminar flow, T restricted to Og satisfies (Gl'). Condition (G2') was verified in Lemma 5.5 
while (G3') and (G4') follow from the analysis of section[3]and section|4](in particular, note that since the 
mixed solution curves obey the s calin g |£| < C\A - A*\, for some C > 0, we cannot have Of = 0.) Finally, 
(G5') is a consequence of Lemma 5.4 since S$ = ^"'({0}) fl Og. We are therefore justified in taking f K' &j to 
be the global extension of C[ ., for i = 1, . . . ,N. Alternatives (i)-(iii) follow from Theorem 6. 1 e) applied 
to 7^. " □ 

Ideally, one would like to eliminate all but the first and second alternatives of Theorem 6.4 as this 
would imply that either the curve is unbounded, or the problem degenerates (i.e. loses ellipticity) along 
Ua>o ^/ ■ As we shall see in section|8] put into the Eulerian framework both of these possibilities lead to us 
to conclude that there is a curve of solutions to the original problem ( |l.l| l-( [L6| along which either the fluid 
tends to stagnation at a point, or the horizontal velocity (in the moving frame) approaches -co somewhere 
within the fluid. 

Typically one excludes the undesirable alternative (iii) by means of a nodal property argument. That 
is, using maximum principles one shows that the sign of certain quantities are invariant along the curve, 
and then argue that this is violated if is a closed loop (or, if we are considering simple bifurcation, if C' s 
returns to T) (see, for example, HIOI ). Taking <x = 0, an argument of this kind was successfully carried out 
in the constant density case (cf. |6|), and, to a lesser degree, in the stratified case (cf. [36]). Unfortunately, 
the introduction of surface tension — in particular, the higher derivatives on the boundary that accompanies 
it — seems to destroy the maximum principle structure of the problem. The rather conspicuous dearth of 
literature on unbounded curves of capillary and capillary-gravity waves is a direct consequence of this 
fact. 



7 Uniform regularity 



In this section we consider the first alternatives of Theorem 5.1 and Theorem 6.4 As in [6|, our goal is the 
establishment of uniform bounds in the Holder norm along the continuum C' g , for the simple eigenvalue 
case, and along the curves { < K' S }-_ 1 for the higher multiplicity case. In that paper Constantin and Strauss 
were able to control the C 3+t? -norm along the continuum in terms sup fteC , sup^7? ; , (assuming that Q is 
likewise bounded along the continuum). However, the addition of the curvature term k[K\ to the nonlinear 
boundary condition complicates matters: we must now allow for the possibility that the curvature may 
blowup (in L°°) along the continuum. 

Theorem 7.1 (Uniform Regularity) Let 6 > be given and assume that fi e C 3+a ((0, \po\)). If 
sup sup /i 7 , < co, sup Q < co and inf inf/cf/z] > -co, 

(Q,h)eC' s r (Q,h)eC' s (Q,h)eC' s T 

then sup (aft)EC , ||/ill c3 „ (5) is finite. 

In the proof of Theorem |7.1| the most technically complicated step is obtaining bounds on the second 
derivatives of h along C' s . The L°° estimates will be a consequence of po being constant along the contin- 
uum, while the higher-order terms can be expressed (via the height equation) in terms of the lower-order 
terms in the usual way. To treat the second-order terms, we shall employ a suite of a priori estimates for 
nonlinear elliptic equations with Venttsel boundary conditions due to Luo and Trudinger (cf. [22]). We 
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shall not, however, need these results in their full generality, so instead we streamline them into a sin- 
gle statement. Of course, in order to apply this to our problem, we will again need to use the method of 
freezing the operator Q. 

First note that due to the periodicity in q, the domain R can be considered as embedded on a torus, al- 
lowing us to ignore the seeming lack of smoothness at the corner points (0, 0), (0, po), (2n, 0), and (2n, po). 
We consider a differential operator F = F(x,z,^,r) e C 2 (R x R x R 2 X S,R) and boundary operator 
G - G(x,z,^, r) 6 C 2 (R X R X R 2 X §, R). Here § denotes the space of 2 X 2 real symmetric matrices, D 
denotes the gradient operator and D 2 the Hessian. Both of these are taken as acting on the space of smooth 
real-valued functions on R which are 27r-periodic in the first variable. As usual, we say that F is elliptic 
at (z,g, r) e R X R 2 X § provided that the matrix F r :- [AF/oVy]^,^ is positive definite at that point. 
Moreover, if A\ and A2 denote the minimum and maximum eigenvalue of F r , respectively, then F is said 
to be uniformly elliptic provided that the ratio A2M1 is bounded. 

We now specialize this framework to the problem at hand. Let a, 7 = ay(z,f), b - b(x,z,^) e C 2 (R X 
R X R 2 , R) and u, 7 eS,b = b(z, f) e C 2 (R X R 2 , R) be given. We consider the case where (F, G) are of the 
form 

F(z, £ r) = otjiz, £)r y + b(x, z, f), G(z, f , r) = a u r u + b(z, 0. 

Here ay is assumed to be such that F is uniformly elliptic, in the sense of the previous paragraph. We shall 
take G to be of Venttsel-type, which means the following: (i) (Uniform obliqueness) There exists some 
X > such that, at each point (q, 0), the inward normal derivative • (0, 1) > x > 0, for all (z, e R X R 2 ; 
(ii) (Ellipticity) uy££j > 0, for all £ = (fi, £2) e K 2 sucn tna t f i s tangent to the lower boundary ofR. 

With F and G given as above we wish to consider the quasilinear elliptic (Venttsel) boundary value 
problem 

F(h, Dh, D 2 h) = 0, in R 

GQi, Dh, D 2 h) = on p = (7.1) 
h = on p = po. 

Theorem 7.2 (Luo-Trudinger) Let h e C 5 (R) be a solution to ( |7.1| l that is In-periodic in the first variable 
and for some constant K > satisfies \h\ + \Dh\ < K in R. Suppose further that for some positive constant 
M the functions F = F(z,£, r) and G — G(z,%) satisfy the following structural conditions: 



A 2 < MA U (7.2) 

m in-i, mi < mama + d, (7.3) 

F rr < 0, (7.4) 

|G|, \G Z \, \G f \, \G ZZ \, \G Z? \, \G & \ < M X , (7.5) 

\d 2 a u \, (1 + |0)13 t 3 f flyl, (1 + \{\ 2 )\d 2 a u \ < MA\ (7.6) 

l&^l, 1^1,1^1(1 + 1^1) < + |^| 2 ) (7.7) 

16 J, \b zx l l*«l, (1 + £1)1**1, d + IfDIM- d + l^l 2 )l%l * MA\{\ + |£| 2 ) (7.8) 

(i + If DIM. + lfDI%l < M*(i + If I) (7.9) 

|D 3 a y |, |£> 3 6| < AfAi (7.10) 

\D\<M X (7.11) 



/or a// (z, £,,r) e R X R 2 X § such that \z\ + If < K. Then there are positive constants p = p(M) and 
C - C(K, M) such that h e Cp^(R) and 
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This is a (weaker) restatement of the results presented in [22]. Also, for consistency with presentation 
of similar theorems earlier in [6| and [36], we have slightly altered notation. We now give a brief sketch 
of the proof. 



Proof of Theorem 7.2 It was already shown in [361 that ( |7.2[ l-( |73| ) together amount to the "natural struc- 
tural conditions" (L-T F1)-(L-T F5) presented in [22 1. We will now show that, unsurprisingly, these further 
imply structure conditions (Luo-T A), (Luo-T Al), (Luo-T A2), (Luo-T B), (Luo-T Bl) and (Luo-T B2) 
of G2J- 

First consider the conditions on F. That (L-T Fl) is equivalent to (Luo-T A) is obvious. Since F(; -,-,r = 
0) = £>(•, •, •), (LT F2) immediately leads to (Luo-T Al). Similarly, if we note that D r F - a,;, it can be 
readily verified that condition (L-T F4) implies the first line of inequalities in (Luo-T A2). To get the 
second line, we appeal to (L-T F3), noting once more that F evaluated at r — is b. 

Next we turn to G. Clearly, (Luo-T B) follows from uniform obliqueness. Since G{-,-,-,r — 0) = 
/?(•, •, •), evaluating (L-T G2)-(L-T G3) at r — yields the growth conditions on /3 in (Luo-T Bl)-(Luo-T 
B2). 

Combining this with the new structural conditions ( |7.6| )-( f77TT| >, the hypothesis of our theorem encom- 
passes all the conditions of 0221 . To arrive at the statement above, we simply apply in succession each of 
the results in §2-§4 of that paper. □ 

Our tools fully assembled, we are ready to prove Theorem UA\ 



Proof of Theorem 7.1 We begin by proving the L°°(7?) norm of h is uniformly bounded along the contin- 



uum. To do so we note that by ( |2.9| l 

< h(q, p)=y + d< tj(0) + d, 
and thus the definition of p$ in ( |1.10| l yields 

17(0) + d < ^ r = \ PX >\ SUp h p 

Thus assuming the hypothesis of the theorem, that is that h p is uniformly bounded in L°° along the contin- 
uum, then the above inequalities together give the boundedness of the L°°-norm of h. 

To bound h q we proceed as in Theorem 6.1 of ||36ll to show that h„ cannot attain an interior minimum 
or maximum. Indeed, since the interior equation here is identical to the one considered in that paper, a 
verbatim application of their argument suffices. As h q = on the bottom, we need only concern ourselves 
with the top. The nonlinear boundary condition then gives 

inf h 2 p (Q - 2cr K [h] - 2gph) < 1 + h\ < sup h 2 p (Q - 2cr K [h] - 2gph). 



T 



The quantity on the left-hand side is bounded uniformly from below by 8 3 on C' 5 . Thus sup^ h q - sup r h q 
can be controlled by sup c , sap-gh p and sup c , sup T (-K[h]). 



For the second-order terms we need to make use Theorem 7.2 and so we shall once again freeze the 
operator Q in the sense of section [5] Fix K > 0, r e R. Let £ = £ (£) be a smooth cutoff function with 
< ( < 1, C = 1 on B K (0) n : £, > S] and supp f C B 2K (0) D : £; > 8/2}. For each x = (x u x 2 ) e R 2 , 
zeR,^ = (^,^ 2 )eR 2 ,put 



b(x,z,Z) := g (/3{-x 2 ) - g{z - T) Pp (x 2 )) 
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and 



Q - 2gp(0)z 



e 2 , 



m + (i - ami- 



Suppose that for some t, (h, Q) is a solution to \1 .\\ with F and G as in Theorem |7.2| that is, 

F{x, z, g r) = (*§n i + (1 + £> 22 - 2f ^ 2 r 12 + (\ (/?(-x 2 ) - g(z - t) Pp (x 2 ))) ((g) 
+m(l-((g)) 



G(z,(,r) = 2crr n + 



2x3/2 



Q - 2gpi0)z ■ 



e 2 



((g+( z (i -((ox 



It follows that each (h, Q) e C' 5 with \h\ + \Dh\ < K satisfies the Venttsel boundary value problem ( |7.1| l for 
the above choice of F and G. The next task is verifying the structural conditions. It is easy to see that F is 
uniformly elliptic. To show uniform obliqueness note that, at any point x e dR, the inward normal is given 
by: 



bf • (0, 1) 



2(1 + 



e 2 



(ig) - (d h ((g))g + 1 - ((g)- 



Inside the region where ( = 1, this is clearly uniformly positive. Also, since we may take ( to be mono- 
tonic decreasing in the ^-direction, the last three terms guarantee that the innward normal derivative is 
uniformly bounded away from zero in the entirety of R 2 . 

The final structural conditions are easy. Since a,;, b, a, ; and b are all suitably smooth, by compactness 
we can choose a sufficiently large M so that ( f7T2]i-(|7.1 \\ hold for ( e B K {$) n {( : ( 2 > 6). (Note that here 
we make use of the assumption y e C 3 to get ( 7.10| >.) 

The more serious step is to gua rantee that h has enough regularity to begin with. That is, we know 
h e C 3+a (R), however Theorem 



7.2 



requires it to be of class C . In order to conclude that C' s lies in C (R) 
we shall apply standard linear elliptic theory. Setting aside the question of f3 for the moment, the fact 
that F is quasilinear tells us that, if h e C k+a (R) with k > 2, then the coefficients of the problem are at 
worst C k ~ l+a (R). Thus, if we assume e C 3+a ((0, \po\)), we have b e C 2+a , which implies by the 
estimates of Lemma 5.2 for the linear Venttsel problem that h e C 4+a (R). This means, in turn, t hat th e 
coefficients are actually C 3+ ". Turning the crank once more, we conclude h e C 5+a (R). Theorem 
therefore applicable, and thus the second order derivatives are bounded. 



7.2 



Thus (F, G) satisfy the structural requirements of Theorem 7.2 We conclude that there exists a con- 
stants C = C(K, M, t) and /j. = p.(K, M, r) so that, for any solution (h, Q) of ( fTTT) with \h\ + \Dh\ < K, 



PI 



< C. 



Note that every h e C' s is such a solution for t = d(h). In order to complete the estimates of the second- 
order terms, therefore, requires removing the dependency on t from C. Toward that end we recall that 
r = d(h) < sup R \h\. Thus we may select K independently of r. Moreover, since in the definition of F, r 
does not occur in any of the second-order terms a,j, it does not affect A\ and A 2 . Indeed, the characteristic 
polynomial satisfied by the eigenvalues is 



A\ - (1 + (\ + %)Ai - (1 + (])( 2 2 ~ *m = 0, ''=1.2 



t-li-2 



which does not depend on r. This is an obvious consequence of the fact that the effects of density variation 
in the height equation are consigned to the lower-order terms. The term in F that does contain r is b, and 
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it can be easily estimated in terms of K. Thus C and p can be chosen independently of r, which gives the 
second-derivative bounds. 

Finally we establish third derivative bounds. Denoting 9 := h q , we differentiate the height equation in 
q to find, 

'(\ + h 2 q )e pp -2h p h q e pq + h 2 p e qq 



C(h, hp 7 Hq, hpp 9 hpq, hqq) 



in R, 



2crh 2 p (0,(1 + h\r^) q + 2h q 6 q + (2gph + 2<TK[h] - Q)h p 6 p 

- D(h q , h p ) on T, 



(7.12) 



9 = 



on B. 



Thus 8 solves a second-order elliptic PDE with Cp^(R) coefficients and right-hand sides C e Cp er (R), 

and D e C^ a {T). Since we have already proved h is bounded uniformly in C^(R), it follows that the 

right-hand side of (J77T2) is bounded uniformly in Cp el (R) along the continuum. As we have seen in section 
[5j so long as h p is bounded uniformly away from along the continuum, problem ( |7.12[ l above will be 
uniformly elliptic with a nondegenerate Ve nttse l-type boundary condition. We may therefore apply the 



linear theory for such problems (e.g. Lemma 5.2k to obtain uniform a priori estimates of 9 in ci^(R) over 



all ofC' s . 

It remains only to bound h ppp in Cp er (7?). By means of the height equation, we may express h pp in 
terms of the lower order derivatives of h: 



h pp = -(l + h 2 q ) 1 {hpy{-p) + h qq h 2 p - 2h q h p h pq ) 



in 



But the right-hand side above is in Cp^(R) by the arguments of the previous paragraphs. Thus h ppp 
Cp er (R) as desired. 

To transition back to the original Holder exponent, a, we merely note that if h e C^(R), then it is 
certainly in C 2 ^(R). The arguments we have used to derive the third derivative bounds in no way relied 
on the particular value of p, so running through them with a instead we find h e C^(R). □ 



8 Proof of the main results 



With the uniform regularity results we established in the previous section, we are now prepared to begin 
the task of proving the main theorems. Theorem 5.1 gave us three possibilities for the continuation of the 



local bifurcation curve(s). In order to prove Theorem 1.2 we need merely to make a careful study of these 



alternatives in light of Theorem 7.1 



Proof of Theorem 1.2 Fix 8 > 0. If C' s is unbounded in R X X, then at least one of the following must 



(i) there exists a sequence (Q n , h n ) e C' s with lim, MOO Q n = oo; or 

(ii) there exists a sequence (Q n , h„) e C' s with lim,,^^ min^ K[h n ] = -oo; or 

(iii) there exists a sequence (Q n , h n ) e C' s with lim,,^^ max-^d p h„ = oo, 

where the second and third possibilities are a consequence of Theqrem |7.1| If C' s contains a point of dOg, 
then either 



(iv) there exists a (Q, h) e C' s with d p h = 8 somewhere in R, or 
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(v) there exists a (Q, h) e C' 5 with 2gph + 2<tk[K\ = Q — 6 somewhere on T. 

Thus, Theorem |5.1| tell us that one of the five alternatives above occurs, or else the bifurcation curve 
intersects T at a point other than (Q*, H*). 

Consider the first alternative (i). By the definition of po in ( |1.10| l, we can estimate 



inf yjp(x,y)(c - u(x,yj) < — - < sup ^p(x,y)(c - u(x,y)\ 
and thus, appealing to (|2.5|), we deduce 



Q = p(0, /7(0))(c - «(0, T](0))f + 2gp(0, T](P))(.ri(0) + d) 

n2 t 2gp(0,77(0))p (8.1) 
<supp(c- M ) + — — -. 

D, 1Ilf -D n <P {C ~ ") 

If for some 6 > the first alternative holds, we conclude from \%.\) that for the corresponding sequence 
(u n , v n ,p n , J] n ) of solutions to ( |1.7| l-( fL8l l, either sup^— ^[p^ic — u n ) — » oo, or inf^ s[fhi(c - u n ) — > 0. 

Similarly, if the third alternative (iii) holds, then, simply by rearranging the change of variables equa- 
tions in ( |2. 1 1 1 >, we have 

d p h n 



Jp~ n (c - u„) 

whence max-^d p h n — > oo if and only if inf^— s/fh&c — u„) — > 0. 

Now consider the second alternative (ii) wherein, along some sequence, we have that a point of infinite 
(negative) curvature is developing. Then, by the dynamic boundary condition ( |1.4| l, we see that 

„ / , ,\-3/2 

Pn ~ Pztm = -OT% (1 + (lj'Xj 

and so max^ n (P„ - P atm ) — > 00 • In other words, as the free surface develops a corner point, the pressure 
at that points (relative to the atmospheric pressure) goes to positive infinity. This, we anticipate, leads to a 
blow-up in the energy on the free surface. Indeed, from the definition of Q in the Eulerian framework, we 
obtain 

Qn =2\E n \ rh - P atm + gp„\t,„d n ) 

= 2(P„- P atm ) + [p„ ((m„ - c) 2 + v 2 ) + 2gp„ (tj„ + <4)J 

> 2 (P n — P a tm) • 



In 



Thus alternative (ii) reduces to alternative (i). 

Next, suppose that for some decreasing sequence 8„ — » the fourth (iv) alternative holds. Again, 
appealing to ( |2.1 1[ ), we see immediately that the corresponding sequence in the Eulerian framework, 
(u„, v„,p„, T] n ), must satisfy sup^j— y[p^{c - u„) — > oo as n — » oo. 

Finally, if we have a sequence 8 n — » for which the fifth alternative (v) holds, then there exists a 
sequence {(Q„,h„)} with (Q„,h„) 6 for each « > 1 and sup r (<2„ - 2gph n 2crK[h n \) — > 0. If we evaluate 
the boundary condition on T for this sequence, we find 

/ 1 \ (l+(d c ,h n ) 2 \ 
max TL , ,o - max / Q , ^9 = max (2« - 2^p„/i„ - 2crK[h n ]) . 



T \(d p h n ) 2 ) T \ (d p h n ) 2 ) T 
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We infer, therefore, that max^<5 ; ,/?„ — » oo, and hence the fifth alternative implies the third. 

By construction, the family of continua C' s indexed by 6 > is increasing as 8 decreases. We may 
therefore define C '.= \Js>qC' s to be the maximal continuum. Then by the considerations of the §2, and 
in particular Lemma 2.1 there exists a connected set C of solutions to ( | 1 . 1 | >-( p~6] > corresponding to C . 
The arguments of the preceding paragraphs, moreover, imply that along some sequence in C, either (i) 



u n — > c, or (ii) infr 



-oo, or (iii) C contains more than one laminar flow solution. This proves 



sup A... 
part (a). 

Part (b) follows from part (a), taking into account the remarks following Theorem 6.4 That is, for each 
6 > 0, we have a path-connected set TCj c C' s where, along some path one of alternative (i)-(v) occurs or 
'Kg is a closed loop. Arguing as before, it is clear that, if 'K' := U<5>o c C' , then there is a corresponding 
path-connected set "K c C of solutions to ( |1.1| |-( |1.6) . It follows that either 7C is a closed loop, or along 
some path in 7C we have either (i) sup^- u —> c, or (ii) infg- u — » -oo. 

Finally, (c) is an automatic consequence of the definition of the space X. The proof is complete. □ 



Proof of Theorem 1.6 Fix £ > and e {1, 2, 3,4}. In view of Theorem 6.4 and Theorem 7.1 the following 
list of alternatives is exhaustive. 

(i) — » oo, as 5 — » oo; or 

(ii) minr /c[f),(i)] — » -oo, as 5 — » oo; or 

(iii) max^d p i)i(s) — * oo, as s — > oo; or 

(iv) min-^d p i)i(s) — > 5, as 5 — » oo; or 

(v) min r (A(j) - 2gph,(s) - 2tr«B,(5)]) 

(vi) "Kgj is a closed loop. 



S, as 5 — > oo; or 



Here (Q,( ), !),(•)) is a global parameterization of 'Ksj in the sense of Theorem 6.1 (by abuse of notation we 
have suppressed the dependency of the parameterization on 6). 

In light of Lemma 12. 11 toTC there is associated a path-connected set 'Ksj of solutions to the problem 
in its Eulerian formulation ( |l.l[ l-( fL"8j ). By the same arguments as in the proof of Theorem 1.2 we deduce 
that the first three alternatives (i)-(iii) imply that one of the following occurs: 



sup u„ 



inf u n 



for some sequence {(u n , v„,p„, r\ n , Q„)} c 'Kgj. The same holds true, again by the considerations of Theo- 
rem 1.2 if either alternative (iv) or alternative (v) hold for a decreasing sequence 6„ with 8 n — > 0. If we 
take "Kj :— IJi5>o ^S,i an d K := IJ; t nen tne theorem is proved. □ 

Remark 8.1 Naively, one might hope that since each < K' B j has an analytic parameterization, the same might 
be true of the associated curves 'Ksj- Unfortunately, this does not seem to be true in general. The difficulty 



lies in Lemma 2. 1 where the key step in transitioning from the height equation formulation back to the 
Eulerian is to reconstruct ijj from (h, Q). In [6,36 1, this is accomplished as follows. Let (Q, I)) be an analytic 
parameterization of 7C^. and, fixing s e R + , let h = l)(s). In view of ( |2.11[ ), for each xq e R, we define 
i//(xq, •; s) to be the unique solution of the first-order ODE: 



^ y (xo,y; s) = -h p l (x , -i^(xo,y, s)) 
iftxo, h(x Q , 0) - d(h); s) = 0. 



(8.2) 



Following the arguments of [6 36], we can show that this ip indeed corresponds to the pseudo-stream 
function for a solution of ( | 1 . 1 1 > — < fT~8^ > , thus obtaining Lemma 2.1 for fixed s. Varying s gives a curve of 
pseudo-stream functions and the path-connected set 'Ksj. But, since t)(s) need not be analytic in p, we 
cannot conclude from (|8.2|l that tfr(-, •; s) is analytic in s. □ 
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9 Examples 

In this section, we consider the implications of our results in a a few of the most important special cases 
for the given functions p and B. 



9.1 Pure capillary waves 

Examining ( |1 .7[ > — ( p~~8] », it is apparent that the contribution of surface tension comes in the form of a term 
involving second-order spatial derivatives, whereas the gravitational constant appears as the coefficient of 
first-order derivative terms. When considering flows in which the wave length is small with respect to <x, 
therefore, we expect the motion to be driven largely by the effects of capillarity. In such circumstances, 
it is often appropriate to discount the presence of gravity entirely, both in the interior of the fluid and on 
the free surface. A common example is that of a pebble being thrown into a quiescent body of water; the 
ripples that result are typically small enough in width as to be accurately modeled as pure capillary waves. 

Strictly speaking, of course, pure capillary waves lie outside of the regime we are considering in this 
paper, as we have taken g to be a strictly positive constant. Nonetheless, our arguments can be easily 
altered to include this scenario. Indeed, defining the pseudostream function and proceeding as before, one 
arrives at a height equation formulation analogous to ( |2.13[ ): 

C (1 + h 2 q )h pp + h qq h 2 p - 2h q h p h pq = -h 3 p B(-p) p Q < p < 0, 

ll+h 2 q +h 2 p (2cn<[h]-Q)=0 p = 0, (9.1) 

[h = P = Po- 



As one would expect — and, indeed, ( |9.1| l confirms — ignoring gravitational effects completely negates 
the significance of stratification. That is, for every solution of ( |9.1| l and any choice of streamline density 
function, there corresponds a solution to (jTTTJ — (JTT8J (cf. 081 1. 

The system of equations ( |9.1| l was first considered by Wahlen in [331, where he used it to give a 
systematic study of small amplitude solutions. In his thesis, Wahlen expanded upon these results, proving 
a global bifurcation theorem in the same vein as Theorem 5.1 (cf. 1351 ). As an easy byproduct of our work 



in sections [6] [7] and [8] we are now able to significantly extend upon this result. 

Let us first begin by presenting the key results for the local theory developed in 113 31 . In keeping with 
our previous notation, we have written these in terms of the Bernoulli function B, but it is worth noting 
that, since p p = for pure capillary waves, B reduces to the familiar vorticity function y of [6,7, 33ll32l 

cm. 

The first step is to establish the existence of a curve of laminar flow solutions, i.e. those where h q = 0. 
In light of the boundary condition on the top, this task is somewhat simpler than for capillary-gravity 
waves. In fact, there is no need to introduce an additional parametrizing variable. 

Lemma 9.1 (Wahlen, 113310 For each Q with < -2B min < Q, the corresponding laminar flow solution 
H = H(p; Q) of fl2T3l is given by 



H{p;Q) := f -^=^==, />„</><(). 

J Pa 



ds 

, n VG + 2fl(5)' 



The existence of small amplitude solutions is proved by a bifurcation argument as was carried out in 
section [3] with the end result being the following. 
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Theorem 9.2 (Wahlen, [33]) Let the wave speed c > 0, the relative mass flux po < 0, coefficient of surface 
tension cr and the Bernoulli function B € C 1+£ *([0, \po\J) and a e (0, 1) be given. Then, if 



sup 

p\£M 



C(2B(p)-2B min )^dp 



> 1, (9.2) 



where 



M : = \ Pl e [p ,0] :crp\> f (p - p x )\2B{p) - 2B min )dp 



>1<\ 

there exists a C 1 -curve T>i oc c of small amplitude traveling wave solutions (u,v,tj, Q) of ( |l.l| l-( fL"6} 
with period 2k, speed c and relative mass flux po , satisfying u < c throughout the fluid. The curve T)\ oc , 
moreover, exhibits the structural properties of Theorem |3 . 1 1 (i)-(iii). 



Note that the statement above differs from the theorem presented in 031 in two ways. First, we have 
specialized to the case where L — 2n. This is done simply as a matter of taste and to stress similarities with 
(6) and 0"6ll . More importantly, we have increased the required regularity of the solutions (u, v, rj) e T)\ oc 
and the vorticity function by one (that this is valid is straightforward to show, by the remarks following the 
statement of \\.2 1.) The motivation for working with more regular solutions is merely to allow us to later 
employ Theorem 



7.1 



Also, let us draw attention to the fact that, for pure capillary waves, the bifurcation from T is always 
simple. This is to be expected as sending g to zero in ( |3.16| l, we find <x c = so that L\ — R + and E% U.Z3 = 0. 
Interestingly, then, double bifurcation is truly a capillary-gravity wave phenomenon: it is absent in both 
pure capillary wave and gravity wave regimes. 

Now, observe that the arguments of section [5] does not in any way rely on assuming g to be strictly 
positive. In fact, by taking g — and replacing the local theory of section [3] with that above, we recreate 
exactly the global bifurcation theorem for pure capillary waves proved by Wahlen in [35]. Let 8 > be 
given and define D' g , Og and Ss as in section |5j but with T>\ oc in place of (?ioc - 

Theorem 9.3 (Wahlen, [35]) Let 5 > be given and let 3Dg be as described above. One of the following 
alternatives must hold: 

(i) D' s is unbounded in Rx X. 

(ii) T)' s intersects T at more than one point. 

(iii) T)' 6 contains a point (Q, h) e dOg. 

Fortunately, our arguments in the remaining sections are similarly indifferent to the strict positivity of ' g. 
For instance, taking g — and running through the proofs of section [7] verbatim yields a theorem identical 



to Theorem 7.1 but for pure capillary waves. One way of interpreting this robustness is to note that, from 
a technical standpoint, the difficulty introduced by stratification and capillarity are neatly confined to the 
interior and boundary, respectively. They must therefore be treated more-or-less independently, and hence 
we have already done the heavy lifting necessary to obtain the corresponding results for the pure capillary 
case. 

Combining these results (and likewise generalizing the proof of Theorem |1.2| i yields the following. 

Theorem 9.4 Fix a wave speed c > 0, wavelength L > 0, relative mass flux po < and coefficient of 
surface tension cr > 0. Fix any a e (0, 1), and let the function B e C 3+a ([0, \po\]) be given such that 
condition ( |9.2| l. Let p e C 1+a ([po, 0]; R + ) be any choice of streamline density function. 

There exists a connected set T> c 5? of solutions (u, v, p, r\, Q) of the traveling, stratified, pure capillary 
wave problem. The set T), moreover, exhibits all of the structural properties described in Theorem [F2](a)— 
(c). 
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This theorem represents a substantial strengthening of Theorem 9.3 elucidating as it does the meaning 
of alternatives (i) and (iii) in the Eulerian framework, as well as proving the existence of a path-connected 
subset. 



9.2 Irrotational homogeneous capillary-gravity waves 

This simplest case to consider with cr > is when f3 = and p = po, a positive constant. Since p p = 0, 
the laminar problem simplifies greatly and we are able to solve it explicitly. In particular, for A > 0, 
po < p < 0, we have 



HM = EZ£>, Q(A) = A + 2gPolPo1 



VI ~ ' \[a 

Thus Aq — (gpo\po\) 2 ^ 3 and a = yA. Furthermore, since we are working in the constant density regime, A c 
must coincide with Aq, whence 

<r c = l(8Po) 2 \Po\ 3 A Q 5 ' 2 = |(^Po) 1/3 IPoI 4/3 - 

Arguing as we have in the previous sections, we see that there exists a local curve of non-laminar 
solutions of period 2n/n bifurcating from T at (//(•; A t ), A t ), for some n e N x , A t > 0, provided there 
exists M = M(p) satisfying 

(a 3 M" - n 2 aM, p Q < p < 0, 
M(p ) = 0, 
Al /2 M'(0) = (gp + nV)M(0). 

From the first equation above we infer that any such M must take the form 

M(p) = sinh| rc (P ^_ P °- > j, p < p < 0. 

Imposing the boundary conditions on the top, it follows that bifurcation will occur for (n,A t ) satisfying 
the following relation 

/U = «V±gPo tanh |jgj (93) 

Note that, for fixed «, there is always a unique choice of A, for which \93\ will hold. This corresponds to 
the fact that, for constant density and zero vorticity, the Capillary-Gravity Local Bifurcation Condition is 
satisfied automatically. However, it may very well be the case that (n\,A„) and (ri2,A t ) each satisfy ( |9.3| l 
for some A„ > and distinct ni,«2 e N x . 

For definiteness, take ti\ — 1 and let A t be chosen so that (1,A») satisfies relation \93\ . If there is no 
«2 > 1 such that ( |9.3| l holds for («2, ^*), then cr e L\ and we are in the simple bifurcation setting described 
by Theorem|3.1| In view of Lemma '. 



. 3.9 we know that simple bifurcation will occur whenever cr > cr c . This 
fact can be seen directly from ( |9.3 1 since, for any such cr, we must have that A is strictly increasing as a 
function of n. 

On the other hand, there is a set E2 of cr lying to the left of cr c for which there will be an n% > 1 such 



setting of Theorem 1.4 As we might expect, the set E2 is highly non-generic. Indeed, the following result 



shows that it has measure zero. 



Proposition 9.5 For homogeneous, irrotational, capillary-gravity waves, E2 is countable. In fact, it con- 
sists entirely of isolated points. 



66 



Samuel Walsh 



Proof Let <x t e E2 be given. Since ( |L-B[ > is always satisfied for irrotational waves, for each cr € R + , there 
exists a unique /I, > eo for which 

-(aVi)' + «pVi = for p Q <p<0, cf>(p ) = 0, ^ /2 0i(°) = + gp(0))«M0). 

has a (nontrivial) solution 0i, where = fii — — 1. By continuous dependence on parameters, we may 
view A t as a smooth function of o~. 

As cr„ is a point of double bifurcation, there exists precisely one other nonpositive integer, P2 — ~ n \, 
for which the above eigenvalue problem has nontrivial solution with p = p2- Clearly ni + 1, since this 
would entail simple bifurcation. Let us first consider the case where «2 ^ 2. 

Again appealing to continuous dependence, we may consider the negative, real eigenvalues p\ and 
jj.2 to be smooth functions of (A, cr) defined in a sufficiently small neighborhood Ki X H2 of (A*,cr*). 
Unravelling definitions, we see that cr e Kj 1S an element of E2 if and only if ^-^(A^cr), cr) is a positive 
integer. It follows that to show Z2 is countable it suffices to confirm 



d(p 2 ) 



* 0. 

(-i„cr») 



do- 

We now prove this claim. Evaluating \93) for n = 1, we find 



A, =(£T + m) )tanh|^Lj. (9.4) 
Differentiating the relation in cr and rearranging terms, we are able to compute 
dA * t u( \ (t , f , s u*( 

^ = tan Mvi:)r (cr+ ^ 0) 2lF sech te 

= A* |cr + gp + (cr + gp ) ^-^sech ^— = 



(9.5) 



which is strictly positive. (In fact, we already knew this in the general case by Lemma 3.10 ) 

Let us denote n(cr) := ^/-// 2 (^*(cr), cr), for each cr e ^ (in particular, observe that n{cr t ) — «2-) From 
( |9.3[ ), therefore, we have 

«/U = (nV + gp ) tanh {^$\ (9.6) 

Differentiating \9.6) by cr, we obtain 



c?n ctt* / dn 2 \ /n|pol\ 
— A* + n—— = \2n— a + n z tanh —7= 

rfcr rfcr \ da \ / 



+ (no- + gp )sech 1 1 ' 



yfTj\yfZdo- 2A 3 J 2 do- 



Upon combining like terms, this simplifies to 
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where 



it. ul n \P®\\ t 2 , \ \Po\ u2 ( n \P0 

C\ :- A t - 2na-tanh| — — - (n cr + gpo)— — sech — — 



dA t 
dcr 



. 2 , ,n\p \ 2 n\p \\dA t 



A J dcr 



All that remains is to prove C2 £ 0. 
Using (|9.5[), we may write 



C 2 = « 



ntanh 



«|/? |\ dA t 



dcr 



1 + (n <T + gpo) TT ^sech 



2A 



— ntanh ( ^= 1 
Returning to ( |9.4[ l and ( |9.6| l, we see 



ntanh |^j l +(n2£r + gpo) j^_ sech 2 



»l/>ol 



tanh 



1 + (o- + gp )-^7rsectf 



jgol 



[C 3 - C 4 ] . 



(ex + gp )sech' 



2/ LpoI 



n/l, sechf ^ )csch( |, 



Thus, C4 can be estimated 



C 4 = 



sinh 



cosh 



Ipol 



sinh 



+ n- 



\Po\ ) 


cosh 


'\Po\ 








"Ipof 


cosh 


'n\po\ 









< 1+ C 3 sinh 2 |MW^ IP01 



Observe, however, that the map 



n t-> sinh 2 ( -^== |csch 2 ( 



(9.7) 



(9.8) 



(9.9) 



has a global upper bound of 1/4 for n > 2 and any values of \po\, A* > 0. Also, since tanh (n\po\/ V^7) > 
tanhflpol/ V^), we have that C 3 > n > 2. From (|9T9j it then follows that C 4 < 3C 3 /4, whence $1$ 
implies C2 > 0. Altogether, then, we have proved that dn/dcr < 0. It follows that the points cr e Z 2 where 
«2 > 2 are isolated. □ 
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Let us investigate further the case when <x e For simplicity, let us suppose «2 > 2 and fix A = A*. 
Note that, by evaluating the right-hand side of ( |9.3| l for n — 1 and « = 2, we find that nondeneracy 
condition (|4.1|l is equivalent to 



4<T + gP0 tanh /2M cotanh / N 



2<x + 2gp \ V^/ 
Now, recalling our earlier notation, we compute, for i — 1,2, 

i "?IpqI . , 2 /»«bol 

- ¥ii = p- - gpo/* sinh — — 

2 VI \ VI 

n?|/?ol Hi 



* 1. 



2VI 2 

j vi r 7jgoi vi 

4 ( 4 4tii 
3n? 



sinh 



/ 2«/bo| \ 

I VI I' 

/ 2»j|pol \ 5VI sjnh / 4n,lPoi y 

I vi / l6n < I vi /, 



+ ^sinh fc^ cosh 3 ' nW 



4 VI 



VI 



VI 



^^4Vl(i^>o|-^VIsinh(^)) 



+ n,n 3 VI 



. / 2(«; - «;)|/?ol 

sinh : 

I VI 

8(«j - «,) 



sinh 



/ 2(fH + Wj)|po| V 
I VI j 

8(n; + «/) 



• u/ n '^°l\ ,l{ n j\Pv\ 

sinh cosh cosh 

2VI \ VI / \ VI / I VI 



(9.10) 



(9.11) 



(9.12) 



Observe that the quantity in ( |9.10| l is strictly negative, as expected. The signs and relative sizes of the 6>,,„ 
and Ojijj are more difficult to ascertain in general, but inserting ( |9.10[ H |9\l2| l into ( |4.2| l yields an explicitly 
verifiable condition under which Theorem [L4]holds. We omit the details. 

To obtain the dispersion relation for these solutions, let U* denote the horizontal velocity corresponding 
to the laminar flow H t . Then, by the definition of po, 

Ipol = J <fc(c-U.)dy = J (d p H.) dy = d^A*. 
Writing \93\ in terms of the depth and the horizontal velocity, therefore, we obtain 



c-V,= 



n z cr + gp Q 



tanh (nd). 
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